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ABSTRACT 


This  manual  is  a guide  for  using-  a computer  program  for  predic- 
tion of  noise  from  freely’  flowing  road  traffic.  The  program 
is  written  in  Fortran  IV  for  use  on  the  DDP-516/H-832  computer 
complex  under  the  control  of  the  Computer  Technology  Division 
of  the  Technology  Directorate  at  TSC. 

The  manual  consists  of  four  parts.  In  the  introduction,  we 
describe  the  limitations  of  the  computer  model.  In  the  second 
part,  acoustical  properties  of  the  traffic  and  sound  propaga- 
tion model  developed  are  presented  and  the  analytical  expres- 
sions used  in  the  computer  program  are  described.  The  third 
part  contains  a description  of  the  structure  of  the  computer 
program  and  of  detailed  calculation  procedures.  The  fourth 
part  gives  practical  guidelines  for  use  of  the  computer  program. 
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1.  INTRODUCTION 


The  relationship  between  road  traffic  characteristics  and  the 
effective  noise  level  at  a given  receiver  location  is  extremely 
complex.  The  characteristics  of  the  sound  emanating  from  a 
particular  roadway  depend  on  parameters  related  to  the  vehicles 
on  the  road,  such  as  vehicle  type,  density,  speed,  acceleration, 
and  directivity  of  radiation.  The  characteristics  of  this 
sound  as  perceived  at  the  receiver  are  affected  by  topological 
considerations  such  as  location,  orientation,  shape,  and  re- 
flectivity characteristics  of  obstacles;  by  atmospheric  and 
ground  absorption;  by  micrometeorological  considerations;  by 
the  formation  of  shadow  zones;  by  interference;  and  by  other 
effects  specific  to  the  particular  receiver  location.  Finally, 
the  predicted  annoyance  level  depends  on  the  rating  scheme 
applied  to  fluctuations  in  the  sound  intensity.  Accordingly, 
complex  models  for  noise  sources,  sound  propagation  paths, 
and  human  psychology  are  required  if  accurate  predictions  are 
to  be  obtained. 


Accounting  for  as  many  of  these  factors  as  possible  should  be 
the  ultimate  goal  in  modelling  road  traffic  noise.  Neverthe- 
less, in  order  to  keep  the  problem  tractable  without  unduly 
sacrificing  predictive  accuracy,  the  following  major  limita- 
tions have  been  incorporated  in  the  current  implementation  of 
the  computer  model; 

a.  The  computer  model  handles  only  freely  flowing  road 
traffic.  Stop-and-go  traffic  with  a time-dependent 
number  of  vehicles  per  length  of  roadway  and  effects 
of  acceleration  on  the  radiated  sound  power  are  not 
included. 


b.  The  computer  model  assumes  a uniform  atmosphere.  Effects 
of  refraction  and  shadow  zones  are  neglected.  Since  the 
sound  pressure  level  at  a distance  of  more  than  a thousand 
feet  from  a source  can  be  strongly  affected  by  micro- 
meteorological  conditions  of  the  atmosphere,  a considerable 
spread  of  actual  data  must  be  expected  around  the  predicted 
values  at  distances  of  more  than  a thousand • feet . 

c.  The  computer  model  regards  all  the  ’contributions  to  the 
sound  intensity  at  a receiver  as  incoherent  and  adds  the 
intensities  without  consideration  of  possible  phase  rela- 
tions. This  simplification  can  result  in  major  errors 
for  the  sound  pressure  level  close  to  absorbing  or  re- 
flecting surfaces. 

d.  In  its  present  form,  the  computer  program  calculates  only 
the  mean  energy  level,  which  is  the  quantity  that  would  be 
measured  by  a special  sound  level  meter  with  a very  slow 
response.  The  calculation  of  noise  data  corresponding  to 
different  rating  scales  for  fluctuating  noise  is  subject 
to  future  improvements  of  the  program. 

Secondary  limitations  include: 

a.  Simplified  assumptions  for  the  orientation  and  shape  of  , 
obstacles  in  the  sound  propagation  path, 

b.  first-order  approximations  • for  the  ’effect  of  ground  cover 
on  the  attenuation  of  sound,  and 

c.  consideration  of  single  reflections  only. 


Despite  these  limitations,  the  computer  program  has  yielded 
predictions  of  A-weighted  noise  levels  that  are  generally 
within  ±3  dB  of  those  measured  in  the  field.  Thus,  the 
approach  used  in  deriving  the  mathematical  models  for  high- 
way noise  appears  to  be  valid,  and  one  might  justifiably 
expect  reasonably  accurate  predictions  over  a range  of 
highway  and  topographical  conditions. 


2.  MODEL  OF  NOISE  FROM  FREELY  FLOWING  ROAD  TRAFFIC 
2 . 1 Energy  Mean  Level 

To  give  the  user  an  overview  of  the  calculations,  we  derive 
briefly  the  mathematical  formulation  for  the  energy  mean  level. 


We  first  consider  the  sound  emanating  from  a single  vehicle 
located  at  a point  x on  a road  segment  extending  from  x^  to 
x^.  The  sound  intensity  1 at  a receiver  located  a distance 
d from  the  road  (see  Fig.  1)  is  given  as 
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^ref^o  ^q-D/10 
xHd^ 


(1) 


where  1 ^ is  the  reference  value  of  the  intensity  at  a short 

distance  r^  (typically  50  ft)  from  a single  vehicle  on  an 
unobstructed  roadway  site,  and  D is  the  uniform  attenuation 
in  dB  of  sound  from  all  points. of  the  road  segment  considered. 
Equation  (1)  is  based  on  the  assumption  of  hemispherical 
radiation  from  a point  source,  for  which  the  intensity  in 
the  diverging  field  decreases  with  the  inverse  square  of  the 
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ROAD  SEGMENT  ^2 


FIG.l  DEFINITION  OF  PARAMETERS  DESCRIBING  A 
ROAD  SEGMENT  AND  A RECEIVER  LOCATION 


fl 


1 I 


I 


] 

] 


Let  us  now  consider  the  situation  in  which  the  road  segment 
x^,x^  contains  a stationary  flow  of  vehicles  having  identical 
noise  characteristics  and  moving  at  essentially  the  same  speed. 
The  mean  vehicle  concentration  A in  vehlcles/f oot , is  given 
as 


A 


(vehicles/ 

V*528o 


(2) 


where  V is  vehicle  speed  in  miles/hour.  The  mean  sound  inten- 
sity at  the  receiver  due  to  all  the  vehicles  on  the  road  seg- 
ment is  then 


I = I 


ref  ° 


10 


-D/10 


A 


xHd' 


dx 


We  define  the  angle  at  the  observer  (Fig.  1)  as 


(3) 


a 


tan  ^ 


X 

d 
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so  that  Eq.  (3)  is  evaluated  as 
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^ref^o  d 


2 A , _-D/10 

3-  Aa  10 


(4) 


where  Aa  is  the  angle  subtended  at  the  receiver  by  the  road 
end  points  x^  and  x^. 

For  T different  groups  of  vehicles,  with  each  group  having  an 
intensity  ^ and  an  average  vehicle  concentration  A^, 

we  obtain  the  mean  intensity 
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AalO 
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t refjt 
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Field  experience  has  shown  that  the  reference  level  L ^ ^ , 

ref , t ’ 

defined  as  ^ = 10  log  ^ dB , has  approximately  a 

normal  distribution  for  certain  large  groups  of  vehicles 
(such  as  passenger  cars  or  heavy  trucks  travelling  at  roughly 
the  same  speed).  Therefore,  rather  than  considering  a large 
number  T of  small  groups  for  all  kinds  of  different  vehicles 
or  for  slightly  different  speeds,  we  consider  a smaller  number 
S of  broader  groupings.  Accordingly,  we  replace  the  summation 
shown  in  Eq.  (5)  as  follows: 


T L „ , /lO 

I 10 

t=l  ^ 


1 
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'L-L 


ref , s 
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ref , s 


10^"^^°  dL, 


(6) 


where  A is  the  average  concentration  of  all  vehicles  of  group 
s,  L „ is  the  mean,  and  a „ is  the  standard  deviation  of 
a normal  distribution  of  reference  levels.  The  expression  for 
the  mean  intensity  from  traffic  with  S groups  of  vehicles  is 
then 


I = — r Aa  10 


-D/10 


S 

I 

s = l 


L . /lO 
10  ref,s 


^(a  . 74.35)^ 

2 ref,s 


(7) 


In  general,  multiple  road  segments  of  a complex  highway  con- 
figuration differ  in  the  distance  d from  the  receiver  perpen- 
dicular to  the  road,  in  the  angle  Aa  enclosed  at  the  receiver 
by  two  lines  towards  the  ends  of  the  road  segments,  and  in 
attenuation  D of  sound  from  the  road  segments.  Thus,  the 
mean  intensity  of  sound  from  complex  road  traffic  is  given 
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by  the  following  double  summation: 
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E 


-D/10 
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segments 


vehicle 

groups 


(8) 


The  energy  mean  level,  in  dB,  is  then  defined  as 


Lg  = 10  log  I 


(9) 


In  the  remainder  of  this  section,  we  describe  the  calculation 
of  the  quantities  employed  in  Eq.  (8)  and  discuss  the  data 
used  in  the  computer  program  for  the  numerical  evaluation 
of  the  energy  mean  level. 

2.2  Reference  Levels 

The  parameters  L ^ and  a ^ defining  the  normal  dlstrl- 
bution  of  levels  at  a distance  r^  from  single  vehicles  of 
group  s generally  depend  on  the  velocity  of  the  vehicles 
and  on  the  frequency  band  considered. 

There  is  no  theoretical  limit  to  the  number  of  vehicle  types' 
that  may  be  considered  by  the  program.  However,  the  current 
implementation  provides  data  blocks  for'  only  two  types  : 
cars  (Type  1)  and  trucks  (Type  2).  ' 

Reference  levels  for  cars  (Lpg^*  stored  in  the  computer 

program  and  are  shown  in  Fig.  2.  A-weighted‘  octave  band  levels 
as  well  as  A-welghted  overall  levels  are  given  for  vehicle 
speeds  of  30  and  70  mph.  The  data  are  based  on  measurements 


(dVa7y^_0L‘Z  3U  GP)  13A31  QNVG 


3AV130  a31H013M-V 
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made  by  N.  Olson  [i].  They  show  the  predominant  influence  of 
tire  noise  with  maximum  levels  at  1000  Hz  for  30  mph  and  at 
2000  Hz  for  70  mph.  An  increase  of  the  A-weighted  overall 
level  by  Ih  dB  for  an  increase  of  the  velocity  V from  30  to 
70  mph  corresponds  closely  with  the  frequently  observed  speed 
dependence  of  40  log  V.  For  velocities  different  from  the 
stored  data,  the  computer  program  interpolates,  using  the 
formula 


L „ , (V)  = L 
ref,l  ^ 


ref  ,1 


(30  mph)  + 


10  log 


_V_ 

30 


10  log 


„„  [L  „ t(70  mph ) 
70  ref ,1  ^ 

30 


- mph)] 


(10) 


Below  1000  Hz,  the  octave  band  levels  at  50  mph  are  about  the 
same  as  those  found  by  Galloway  et  al  [2].  Above  1000  Hz, 
Galloway's  data  are  lower  by  5 to  10  dB. 

In  contrast  to  our  speed  dependence  with  38  log  V for  the  A- 
weighted  overall  level,  Galloway  found  only  30  log  V.  However, 
the  differences  are  only  1.2  dB  at  70  mph  and  1.8  dB  at  30  mph. 

The  standard  deviation  a „ . was  based  on  Olson's  data. 

ref,l 

Although  there  is  some  dependence  on  frequency  and  velocity, 
we  propose  to  use  the  single  number  -^  = 2.^  dB.  This 

number  is  consistent  with  findings  in  the  USA  [2]  and  in 
England  [S].  Note  that  any  number  smaller  than  3 dB  increases 
the  energy  mean  level  by  less  than  1 dB. 

Reference  levels  for  trucks  (L  „ are  plotted  in  Fig.  3* 

The  data  are  based  on  measurements  made  by  Olson  [I]  and  on 
recent  measurements  at  the  New  Jersey  Turnpike.  We  ignore 
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the  speed  dependence  observed  in  the  higher  frequency  bands  for 
trucks  with  cross-bar  tires,  and  we  account  for  engine,  fan, 
and  exhaust  noise  exclusively  with  speed-independent  levels. 

The  A-weighted  overall  level  is  higher  by  6 dB  than  data 
reported  by  Galloway  [2].  We  found  out  that  Galloway's  data 
are  referred  to  trucks  with  good  mufflers,  but  that  our  data 
are  more  representative  for  actual  truck  noise.  The  influence 
of  mufflers  on  the  A-level  indicates  that  most  of  the  truck 
noise  is  radiated  from  the  exhaust.  The  effective  source 
height  of  trucks  should  therefore  be  chosen  about  8 ft  higher 
than  the  source  height  of  passenger  cars. 

A value  of  3-5  dB  for  the  standard  deviation  a „ „ is  con- 

ref, 2 

sistent  with  data  from  Galloway  [2],  except  for  steep  gra- 
dients of  the  road  where  considerably  higher  numbers,  up  to 
6.8  dB,  have  been  measured. 

2.3  Attenuation 

2.3.1  Atmospheric  absorption 

For  outdoor  sound  propagation,  the  attenuation  by  atmospheric 
absorption  at  a temperature  of  68°F  can  be  reasonably  well 
approximated  by  [4] 

Da  = 7.^  10"®  dB  , . (11) 

where  f = geometric-mean  frequency  of  band,  Hz 
^ = relative  humidity  of  the  air,  % 
r = distance  between  source  and  receiver,  m. 


11 


The  computer  program  uses  a relative  humidity  of  55^  and  octave 
bands  with  the  numbers  n = 2,..., 9 for  mid-frequencies  from 
62.5  Hz  to  8000  Hz.  Distances  are  given  in  feet.  Thus, 

Eq.  11  becomes 


Da  = 7.^ 


.n 


3.28 


::  10 


- 7 


,n 


(12) 


where  r’  is  taken  by  the  computer  program  as  the  minimum  dis- 
tance, in  feet,  between  the  receiver  and  a road  segment.  At 
a frequency  of  500  Hz,  where  n = 5»  the  attenuation  is  less 
than  1 dB  up  to  .distances  of  10,000  ft.  Since  the  attenuation 
of  traffic  noise  is  typically  characterized  by  the  attenuation 
in  the  500~Hz  octave  band,  atmospheric  attenuation  has  hardly 
any  influence  in  all  practical  cases  of  noise  prediction  within 
a distance  of  a few  thousand  feet  from  a road.  Consequently, 
no  checks  are  made  in  the  computer  program  for  small  differ- 
ences in  attenuation  by  atmospheric  absorption  of  sound  from 
different  locations  on  a road  segment;  it  is  assumed  that  the 
attenuation  from  the  nearest  point  on  the  road  is  representa- 
tive for  the  entire  road  segment. 


Note  that  the  attenuation  predicted  by  Eq.  (12)  will  be  too 
low  if  the  actual  humidity  of  the  air  is  considerably  less 
than  55/^  or  if  the  actual  temperature  is  between  20  and  30°F. 
Conversely,  the  predicted  will  be  too  high  for  relative 
humidities  much  greater  than  55^  [4]» 


2.3.2  Diffraction 

Diffraction  of  sound  is  caused  by  obstacles  in  the  direct  or 
reflected  propagation  paths  from  the  roadway  to  the  receiver. 
Such  obstacles  can  be  artificial  barriers,  earth  berms,  hills. 
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buildings,  etc.  For  the  calculation  of  attenuation  by  diffrac- 
tion, we  assume  that  the  obstacle  can  be  modeled  by  a rigid, 
impervious  screen  oriented  perpendicular  to  the  ground  plane 
and  that  sound  is  diffracted  over  the  top  edge  of  the  screen 
exclusively.  We  neglect  the  shape  of  hills  and  the  thickness 
of  barriers  because  of  the  lack  of  available  knowledge.  The 
sound  absorption  and  transmission  properties  of  barriers  are 
not  considered  because  they  play  a minor  role  in  most  practical 
cases.  We  expect  that  the  neglect  of  diffraction  around  the 
ends  of  barriers  v/ill  Introduce  no  significant  errors,  and 
it  simplifies  considerably  the  computational  procedures. 
Furthermore,  a diffracting  barrier  is  now  completely  specified 
by  the  coordinates  of  the  two  end  points  of  the  top  line. 

The  attenuation  of  sound  by  barriers  is  determined  primarily 
by  the  difference  6 between  the  path  length  of  the  shortest 
ray  from  the  source  over  the  top  edge  of  the  barrier  to  the 
receiver  and  the  path  length  of  the  direct  ray  from  the  source 
to  the  receiver  in  the  absence  of  the  barrier  (see  Fig.  ^). 

For  large  path  length  differences,  the  attenuation  in  the 
acoustical  shadow  zone  of  a barrier  is  limited  by  effects  of 
refraction  and  scattering  of  sound  in  the  atmosphere.  Based 
on  data  from  Rathe  [5],  the  computer  program  has  a built-in 
maximum  attenuation  of  24  dB. 

The  attenuation  is  not  zero  for  zero  path  length  difference 
(i.e.,  for  a ray  grazing  over  the  barrier).  For  this  situation, 
the  theory  of  Fresnel  diffraction  yields  an  attenuation  of 
about  5 dB.  The  attenuation  becomes  negligible  when  a direct 
sound  ray  travelling  from  the  source  to  the  receiver  passes 
far  over  the  top  edge  of  the  barrier.  To  simplify  computations. 
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DEFINITION  OF  THE  PATH  LENGTH  DIFlFERENCE  6 FOR 
SOUND  DIFFRACTION  BY  A BARRIER 


we  neglect  diffraction  effects  when  the  height  difference 
between  the  direct  ray  and  the  top  of  the  barrier  is  greater 
than  20  feet. 

For  height  differences  of  20  feet  or  less,  we  calculate  the 
Fresnel  number 


(13) 


where  cS  is  the  path  length  difference  and  A is  the  wavelength 
corresponding  to  the  center  frequency  of  an  octave  band.  We 
assume  normal  atmospheric  conditions,  with  a speed  of  sound  c 
of  1120  ft/sec.  Thus,  for  a center  frequency  f,  the  Fresnel 
number  is  taken  to  be 

N = — 6 . (14) 

c 


We  employ  the  Fresnel  number  and  an  analytic  approximation  to 
measured  data  from  Maekawa  for  the  calculation  of  the  barrier 
attenuation  [4] 


20  log 


/2ttN 


tanh  /2ttN 


+ 5 dB 


for  N > -0.2 


(15) 


otherwise 


The  formula  Eq . (15)  is  applicable  to  both  positive  and  nega- 

tive values  of  N.  However,  for  the  actual  computation  we 
separate: 
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for  N < -0.2 


20  log 


/2tt  I N ~ 
tan/27T  I N 


+ 5 dB 


20  log  + 51  dB 

tanh  /2ttN 


2H  dB 


for  -0.2  < N £ 0 

for  0 < N £ 12.5 
for  N > 'l2. 5 


(16) 


The  last  line  in  Eq.  (16)  accounts  for  the  above-mentioned 
upper  limit  to  barrier  attenuation. 

As  shov;n  in  Fig.  5j  the  attenuation  of  the  A-weighted  sound 
pressure  level  of  typical  passenger  car  noise  (see  Fig.  2) 
is  almost  identical  v/ith  the  sound  attenuation  in  the  500  Hz 
band.  Hence,  the  only  number  important  for  the  attenuation 
of  road  traffic  noise  is  the  path  length  difference. 


The  path  length  difference  accounts  for  heights  and  distances 
of  a point  source,  a receiver,  and  the  top  edge  of  a barrier. 
Furthermore,  it  allows  us  to  account  for  the  reduced  attenuation 
of  rays  oblique  to  the  top  edge  of  the  barrier  [5]. 

For  noise  from  a road  segment  and  for  a barrier  at  oblique 
angle  to  the  road,  the  computer  program  finds  the  path  length 
difference  6,,  for  sound  from  the  nearest  point  on  the  road. 

Then,  by  assuming  a monotonic  variation  of  the  path  length 
difference  from  other  points  on  the  road,  the  extreme  ends 
of  the  road  segment  are  considered.  If  the  path  length 
differences  6 and  6^  for  these  end  points  differ  from  6^ 
by  more  than  a number  that  results  in  an  attenuation  differ- 
ence of  about  1 dB,  the  road  segment  between  the  near  point 
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FIG. 5 BARRIER  ATTENUATION  AS  A FUNCTION  OF  THE  PATH  LENGTH  DIFFERENCE 
FOR  A FREQUENCY  OF  500  Hz  AND  FOR  A-WEIGHTED  LEVEL  OF  TYPICAL 
PASSENGER.  CAR  NOISE 


N and  the  point  1 or  the  point  2,  respectively,  is  cut  in  half. 
New  path  length  differences  are  calculated  for  the  new  end 
points  of  the  road  segment,  and  the  procedure  of  reducing  the 
length  of  the  road  segment  is  repeated  until  the  attenuation 
by  diffraction  is  approximately  constant. 

The  criterion  used  for  the  acceptance  of  a sufficiently  small 
difference  in  path  length  differences  is : 


The  numbers  are  'based  on  a frequency  of  500  Hz,  for  which  the 
effect  of  Eq . (17)  on  the  attenuation  is  plotted  in  Pig.  6. 

In  case  of  multiple  diffraction  by  several  barriers  in  parallel, 
we  account  for  the  strongest  diffraction  exclusively.  This  is 
a conservative  procedure  resulting  in  attenuations  that  are 
somewhat  too  small,  but  it  seems  to  be  the  most  reasonable 
way  to  bypass  the  very  complicated  and  not  yet  fully  under- 
stood problem  of  multiple  diffraction. 

2.3.3  Ground  absorption 

Ground  attenuation  is  a function  of  the  structure  and  the 
covering  of  the  ground,  both  of  which  influence  its  acoustic 
properties,  and  of  the  heights  of  the  source  and  receiver 
above  the  ground. 

In  the  computer  program,  we  consider  in  a very  simple  approxi- 
mation rectangular  ground  strips,  which  are  defined  by  two 
end  points  of  a center  line  and  by  a width,  and  which  have 
either  a low  cover,  such  as  shrubbery  and  thick  grass,  or  a 
high  cover,  such  as  trees. 


(17) 
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FIG. 6 BARRIER  ATTENUATION  AS  A FUNCTION  OF  THE  FRESNEL  NUMBER  OF  DIFRACTIONS. 
RESULTS  OF  MAEKAWA'S  MEASUREMENTS,  • ( 1 6 ) WITH 


The  height  of  a sound  ray  travelling  from  the  source  to  the 
receiver  over  the  ground  strip  is  checked  only  with  respect 
to  the  center  line  of  the  strip.  Thus,  it  is  assumed  that 
the  plane  of  the  ground  strip  is  approximately  parallel  to 
a plane  defined  by  a road  segment  and  a receiver.  If  the 
height  of  the  direct  sound  ray  from  the  source  to  the  re- 
ceiver is  more  than  10  ft  above  a ground  strip  with  a lov; 
cover  or  more  than  30  ft  above  a ground  strip  with  high 
cover,  any  sound  attenuation  due  to  ground  absorption  is 
neglected.  The  heights  of  10  and  30  ft  are  based  on  rough 
estimates  rather  than  on  field  experience  and  might  be 
revised  if  found  necessary. 

In  general,  the  amount  of  ground  attenuation  cannot  be  stated 
in  terms  of  excess  attenuation  per  unit  of  distance.  To  a 
first  approximation,  however,  such  behavior  can  be  assumed 
in  the  range  of  distances  of  200  to  2000  ft  unless  the  total 
attenuation  exceeds  30  dB. 


No  attempt  has  been  made  to  calculate  accurate  distances  over 
a ground  strip  with  the  computer  program.  Instead,  mean  path 
lengths  r over  ground  strips  ahe  calculated  with  the  formula 


r 


tt/2 

i + i 

w I 


(18) 


where  w is  the  v;idth  of  the  strip  and  9J  the  length  of  the 
center  line. 

The  following  analytical  approximations  to  average  values  of 
measured  data  [4]  are  used  to  calculate  the ‘attenuation  of 
sound  propagating 


20 


"k' 


i", 


r^’ 

dihi'- 

yi4: 

[5|i  il4't  ■^y 


5!t#. 


: _ ^ 


:-A- 

■’^0 


?•  -..  "'i. «•' 


mMM^. 


.S3  - X, 


)w  •'■■  ■ 


Wv.:;  ■■:■:■■'; 


§M' 


" -;  ^ ' ffl  . ,.  gj>„  V 


--':lC';''V-V'- 


.■;y=::  ' ♦'fly. .■■;•• ‘.J 

■ ••  ■■''''-  ''Avy 


a.  through  shrubbery  and  over  thick  grass 


Dq  = [0.18  log  (f)  - 0.31]  dB  , (19) 

(see  Fig.  7) 
b.  through  tree  zones 

Dq  = 0.01  dB  ■ (20) 

(see  Fig.  8),  where  r is  given  in  feet  and  f in  Hz. 

For  the,  computer  program,  we  calculate  the  following  attenua- 
tions with  r in  feet  and  with  octave  band  numbers  n,  where 
n = 5 for  the  octave  band  center  frequency  f = 500  Hz: 

a.  for  shrubbery  and  thick  grass 

K0.0l6n  - 0.028)r  dB  if  £ 30  dB 

= 

(30  dB  if  > 30  dB  (21) 

u 


for  tree  zones 


2^/3 

1,310 

30  dB 


r dB 


if  D„  < 30  dB 

U 

if  D„  > 30  dB 

(j 


(22) 


The  notation  in  the  second  lines  of  Eqs . (21)  and  (22)  means 

that  the  attenuation  is  limited  to  30  dB. 

Consistent  with  the  assumption  that  the  ground  attenuation  is 
proportional  to  the  mean  path  length  of  the  sound  over  the 
ground  strip,  the  computer  program  accumulates  attenuations 
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of  various  ground  strips  in  series.  However,  since  the  path 
length  considered  represents  a statistical  average  for  rays 
propagating  in  all  directions  over  the  strip,  the  path  length 
over  two  equal,  parallel  strips  is  not  just  twice  the  path 
length  over  a single  strip  of  twice  the  width  but  is  generally 
shorter.  Consequently,  the  attenuation  calculated  for  one 
strip  will  be  smaller  than  the  attenuation  calculated  for 
two  strips  in  parallel  each  having  half  the  width. 

The  purpose  of  using  the  statistical  formulation  for  path  length 
given  by  Eq.  (18)  is  to  obtain  reasonable  predictions  for  the 
effects  of  ground  absorption  on  the  average.  There  exist, 
hov/ever,  particular  cases  where  the  model  will  not  be  very 
accurate.  For  example,  the  attenuation  of  sound  from  a short 
road  segment  by  long,  narrow  absorbing  strips  is  overestimated, 
whereas  the  attenuation  by  a wide  strip  oriented  perpendicu- 
larly to  the  road  is  underestimated.  To  some  extent,  we 
expect  these  modeling  errors  to  compensate  for  one  another 
in  most  practical  situations.  In  general.  Inaccuracies  are 
inherent  to  the  entire  problem  of  ground  absorption. 

2.4  Reflection 

The  sound  field  at  a receiver  results  from  contributions  of 
direct  (or  diffracted)  and  reflected  rays.  In  many  practical 
cases  of  sound  propagation  from  a highway,  corrections  applied 
for  reflections  are  small  compared  to  the  inaccuracies  involved 
in  the  prediction  of  ground  attenuation  and  in  uncertainties 
with  acoustical  shadow  zones  owing  to  wind  and  temperature 
gradients  in  the  atmosphere.  Therefore,  the  computer  program 
has  been  designed  to  account  for  reflections  with  a first- 
order  approximation. 
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In  the  computer  program,  we  disregard  phase  relations  between 
the  various  contributions  and  consider  incoherent  waves  for 
which  the  total  sound  intensity  is  made  up  by  the  sum  of  the 
intensities  of  the  individual  contributions. 

There  are  always  reflections  from  the  road  surface.  However, 
the  contributions  from  these  reflections  are  already  accounted 
for  in  reference  levels  at  a short  distance  from  individual 
vehicles  on  the  road.  They  don’t  enter  again  into  the  cal- 
culations . 

Reflections  at  the  ground  plane  further  out  from  the  roadway 
are  disregarded  because  they  generally  result  in  a complex 
interference  pattern  with  the  direct  ray.  Consideration  of 
these  effects  is  beyond  the  scope  of  a first-order  approxi- 
mation . 

Reflections  at  any  inclined  plane  result  in  rays  directed 
either  towards  the  ground,  and  thus  being  neglected,  or 
towards  the  sky,  and  thus  not  contributing  to  the  sound 
intensity  at  a normal  receiver  location  close  to  the  ground. 
Therefore,  we  consider  only  reflections  on  planes  that  are 
perpendicular  to  the  ground  plane. 

Within  the  first-order  approximation,  we  further  neglect  the 
actual  frequency-dependent  magnitude  of  reflection  coeffi- 
cients and  distinguish  only  between  reflection  coefficients 
0 and  1 of  reflecting  surfaces. 

In  order  to  determine  whether  a reflective  barrier  is  high 
enough  to  be  effective,  we  consider  the  possibly  reflected' 
ray  that  travels  a minimum  distance  from  the  road  segment 
to  the  receiver.  A reflective  plane  perpendicular  to  the 
ground  is  considered  high  enough  if  the  ray  strikes  the 


barrier  at  least  two  feet  below  the  top  edge  of  the  barrier. 
For  reflection  points  within  two  feet  of  the  top  edge, 
diffraction  effects  are  considered  to  be  strong  enough  for 
all  frequencies  so  that  the  reflected  ray  is  sufficiently 
reduced  in  amplitude  to  be  negligible. 


Also  neglected  are  reflections  from  planes  that  are  either  very 
short  or  very  remote  so  that  the  contribution  to  the  sound 
Intensity  at  the  receiver  is  less  than  10^  of  the  intensity 
received  via  direct  (or  diffracted)  rays  from  the  road  segment 
under  consideration.  The  analytical  formulation  for  this 
criterion  is 


10  ^b/10  < 0.1 

d ' Aa 


(23) 


where  d 


distance  from  the  receiver  to  the  road  segment. 


Aa  = aspect  angle  of  the  road  segment  at  the  receiver, 

d'  = distance  from  the  road  segment  to  a receiver 
location  Imaged  about  the  reflector. 

Aa'  = aspect  angle  of  the  barrier  at  the  image  receiver, 

Dg  = attenuation  in  dB  by  diffraction  of  the  direct  ray 
due  to  a possible  barrier  (referred  to  a frequency 
of  500  Hz) . 


Single  reflections  are  considered  exclusively;  contributions 
from  rays  that  strike  two  or  more  reflectors  are  ignored.  It 
is  essential  for  the  future  calculation  of  higher  order  sta- 
tistical parameters  of  road  traffic  noise  that  the  reflections 
of  sound  from  a certain  road  segment  be  treated  as  amplifica- 
tions of  the  direct  (or  diffracted)  rays  and  not  as  uncorrelated 
contributions  from  independent  road  segments.  The  factor  P 
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multiplying  the  intensity  of  the  direct  sound  from  a road  seg- 
ment is  calculated  using 

F = 10  ^b/IO  ^ Dj^/lO  ^ ^2i|) 

i=l  Audi 
X 

where  the  subscript  i indicates  reflections  at  N different 
surfaces,  each  of  which  might  be  diffracted  by  a barrier  before 
or  after  reflection  and  therefore  might  have  an  attenuation  . 
The  factor  F is  calculated  as  a function  of  frequency.  The 
notation  in  Eq.  (2^)  is  the  same  as  in  Eq.  (23)  except  for 
the  angle  Aa'  which  denotes  the  aspect  angle  of  the  road  seg- 
ment at  the  image  receiver  (see  Fig.  9)- 

2.5  Combination  of  Attenuation  and  Reflection 
2.5.1  Atmospheric  absorption 

The  computer  program  accounts  for  atmospheric  absorption  in 
combination  with  barrier  diffraction,  ground  absorption  and 
reflections.  The  path  length  used  for  calculating  the  atmo- 
spheric absorption  is  the  direct  distance  from  the  source  to 
the  receiver  and  is  not  corrected  for  the  path  length  difference 
6 of  diffracted  rays  nor  for  the  increased  path  length  of  re- 
flected rays.  The  factor  F,  defined  by  Eq . (2^),  is  multi- 
plied by  10  where  is  defined  by  Eq.  (13)  j in  order 

to  calculate  for  each  individual  road  segment  s the  factor 

10“^s/10  = B iq-Da/10  io“^g/10  ^ ^25) 

which  is  employed  for  the  calculation  of  the  energy  mean  level 
in  Eqs . (8)  and  (9 ) • 
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ROAD  SEGMENT 


RECEIVER 


FIG. 9 EXAMPLE  FOR  RAY  TRACES  FROM  A REFLECTOR, 
INDICATING  THE  CONSTRUCTION  BY  MEANS  OF 
AN  IMAGE  RECEIVER 
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2.5.2  Ground  absorption 


The  computer  program  accounts  for  ground  absorption  only  if  the 
rays  travelling  over  the  ground  are  not  diffracted  by  a barrier. 
This  concept  is  based  on  field  experience  with  low  ground 
cover  [2jS].  For  high  ground  cover,  it  might  be  necessary  to 
reconsider  the-  approach  taken.  As  shown  with  the  formula 
Eq.  (25),  the  same  ground  absorption  is  as.sumed  for  the  direct 
rays  and  for  all  reflections. 

2.5.3Reflection 

The  computer  program  accounts  for  reflections  in  combination  . 
with  diffraction  provided  that  there  is  only  a single  diffrac- 
tion before  or  after  the  reflection  and  that  the  path  length 
increase  due  to  diffraction  is  less  than  12.5  ft.  Doubly 
diffracted  reflections  are  neglected  as  well  as  very  weak 
single  reflections  that  suffer,  in  the  500  Hz  band,  the 
maximum  attenuation  of  24  dB  assumed  for  barrier  diffraction. 

The  attenuation  of  reflected  rays  by  diffraction  is  calculated 
for  one  location  on  the  road  segment  only:  the  point  nearest 

to  the  image  receiver.  No  attempt  has  been  made  to  refine 
this  calculation  by  checking  for  the  attenuation  from  other 
points  on  the  road  segment,  since  the  contribution  of  dif- 
fracted reflections  will  be  generally  small  and,  hence, 
inaccuracies  of  the  calculation  will  be  negligible. 
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3.  DESCRIPTION  OF  THE  COMPUTER  PROGRAM 


3 . 1 Introduction 

The  computer  program  is  designed  to  run  on  the  DDP-516/H- 832 
computer  complex  which  is  currently  Installed  in  the  Computer 
Technology  Division  at  TSC.  The  program  is  written  in  FORTRAN-IV 
and  is  designed  to  run  in  batch  mode.  Because  of  the -modular 
organization  of  the  program,  only  a modest  level  of  effort  should 
be  required  to  modify  it  for  an  interactive  mode  of  operation. 

As  currently  implemented,  inputs  are  provided  via  punched  cards 
and  outputs  are  provided  on  the  user  teletype. 

To  meet  the  constraints  on  core  storage  imposed  by  the  PDP-516, 
the  disk  is  used  whenever  feasible  for  storing  large  temporary 
arrays  and  input  data.  Should  the  program  be  implemented  on  a 
computer  which  has  a larger  core  memory,  computation  time  will 
be  significantly  reduced  by  using  core  memory  for  storing  certain 
arrays  nov;  stored  on  the  disk. 

The  organizational  scheme  of  the  program  is  shown  in  the  follow- 
ing flow  diagram.  The  major  blocks  of  code  are  contained  in  the 
main  program  and  in  the  subroutines  INPUT  and  GEOMRY.  The  sub- 
routine INPUT  reads  input  data  provided  on  cards,  stores  this 
information  on  disk  files,  and  outputs  it  on  the  user's  teletype. 
The  subroutine  GEOMRY  provides  the  bulk  of  the  calculations 
performed  by  the  computer  program  and  uses  a number  of  addi- 
tional subroutines  and  functions  in.  the  process.  Additional 
subroutines  called  by  the  main  program  are  FILES,  which  assigns 
device  numbers  to  the  various  I/O  devices  and  files  used  by  the 

i 

program,  and  INTER,  which  computes  reference  intensities  for 
the  various  vehicle  groups  specified  by  the  user. 
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The  main  program  and  the  subroutines  INPUT  and  GEOMRY  are  dis- 
cussed individually  in  some  detail  in  the  remainder  of  this 
section  and  the  remaining  subroutines  and  functions  are  de- 
scribed very  briefly.  Listings  for  the  entire  computer  program 
are  given  in  Sec.  3-6.  Details  on  program  usage  are  provided 
in  Sec.  4. 

The  following  definitions  are  used  throughout  the  remainder  of 
this  manual.  A "road"  or  "roadway"  is  defined  as  a set  of 
connected  straight  "road  sections"  for  which  the  traffic  condi- 
tions are  uniform.  That  is,  the  number  of  vehicle  groups  and 

the  concentration  and  mean  vehicle  speed  associated  with  each 

« 

group  are  the  same  for  all  sections  of  a given  road.  Barriers 
are  treated  similarly.  Thus,  a "barrier"  is  defined  as  a set 
of  connected  straight  "barrier  sections"  for  which  the  absorptive/ 
reflective  properties  are  the  same. 

3.2  Main  Program 
3.2.1  Overview 

The  primary  function  of  the  main  program  is  to  control  the  flow 
of  operations.  Various  subroutines  are  called  upon  to  perform 
the  bulk  of  the  input-output  operations  and  the  analysis.  In 

» 

addition,  certain  variables  are  initialized  in  the  main  program 
and  the  final  computation  of  energy  mean  levels  is  made. 

A flov7  diagram  of  the  main  program  is  given  below.  This  diagram 
is  self-explanatory.  Whenever  possible,  statement  numbers  are 
shown  so  that  reference  may  be  made  to  the  corresponding  block 
of  code  in  the  program  listing. 
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Main  Program 
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3.2.2  List  of  variables  for  main  program 

Variables  appearing  in  the  main  program  are  listed  below.  Array 
variables  are  not  generally  identified  as  such;  the  dimension- 
ality of  the  variables  can  be  determined  from  the  program 
listing  in  Sec.  3-6.  The  variables  XRl  and  XR2  listed  below 
are  vector  quantities  which  indicate  the  coordinates  of  road 
end  points.  The  first  element  of  the  vector  indicates  the 
x-coordinate , the  second  element  the  y-coordinate , and  the 
third  element  the  z-coordlnate . This  format  is  used  widely 
throughout  the  remainder  of  the  entire  computer  program  to 
denote  coordinates. 


CO 

Factor  accounting  for  standard  deviation  of  reference 

level 

I 

Index 

IBAR 

Barrier  number 

IBR 

Barrier  file  device  number 

IDF 

Barrier  subfile  device  number 

lERR 

Error  index 

IGA 

File  device  number  for  absorptive  ground  strips 

IGRA 

Ground  strip  number 

IGTA 

Subfile  device  number  for  absorptive  ground  strips 

II 

Index 

IQ 

Index 

IRD 

Roadway  file  device  number 

IREF 

Reflector  file  device  number 

IRFDF 

Barrier  subfile  device  number 
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ISEG  Barrier  section  number 
J Index 

M Roadway  number 
N Road  section  number 
NB  Number  of  barriers 
NF  Number  of  frequency  bands 
NG  Number  of  absorptive  ground  strips 
NQ  Number  of  vehicle  types 
NQQ  Number  of  groups  within  one  vehicle  type 
NQS  Vector  notation  for  number  of  vehicle  groups 
NR  Number  of  roadways 
NRG  Number  of  receivers 
NRSEC  Number  of  sections  for  one  roadway 
RDIN  Vector  notation  for  initialization  parameters 
SIG  Standard  deviation  of  reference  level 
VEXPH  Vehicles  per  foot 

XLA  A-weighted  overall  intensity  and  level 
XLE  A-weighted  intensity  and  level  in  frequency  bands 
XMPH  Vehicle  speed  in  mph 
XRl  Road  section  initial  point 
XR2  Road  section  end  point 
XRC  x-coordlnate  of  receiver 
YRC  y-coordinate  of  receiver 
ZRC  z-coordlnate  of  receiver 
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3.3  Subroutine  INPUT 


3.3.1  Overview 

This  subroutine  reads  the  input  data  provided  by  the  user  on 
punched  cards,  stores  this  data  on  disk  files  in  a format  suit- 
able for  use  in  the  remainder  of  the  computer  program,  and 
outputs  the  data  on  the  user  teletype  so  that  a permanent  record 
is  created. 

The  computer  program  is  organized  in  such  a way  that  the  user 
need  not  redefine  the  entire  situation  each  time  a new  case  is 
analyzed.  Instead,  it  is  necessary  to  specify  only  those 
aspects  of  the  problem  that  are  different  from  the  previous 
case.  Accordingly,  each  block  of  input  data  is  preceded  by  a 
control  card  that  indicates  the  type  of  data  (e.g.,  barrier 
specifications)  and  the  number  of  items  (e.g.,  number  of  bar- 
riers) to  follow.  After  a special  control  character  has  signi- 
fied the  end  of  data  input,  the  entire  set  of  input  variables 
is  typed  out  and  control  returns  to  the  main  program. 

A flow  diagram  of  this  subroutine  is  given  below,  followed  by 
a list  of  variables. 
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3.3.2 

List  of  variables  for  subroutine  INPUT 

ALPHA 

Alphanumeric  information 

BG 

Width  of  absorptive  ground  strip 

I 

Index 

II 

Number  of  items  in  data  block 

12 

Dummy  variable 

lA 

Alphanumeric  "A" 

IBR 

File  device  number  for  barriers 

IDF 

Subfile  device  number  for  barriers 

IDN 

Index  for  program  initialization  parameter 

IDUM 

Index  for  type  of  absorptive  ground  cover 

IG 

Alphanumeric  "G" 

IGA 

File  device  number  for  absorptive  ground  strips 

IGO 

Index  for  data  blocks 

IGTA 

Subfile  device  number  for  absorptive  ground  strips 

ILAST 

Alphanumeric  indicator  for  last  section 

IR 

Alphanumeric  "R" 

IRD 

File  device  number  for  roadways 

• I REF 

File  device  number  for  reflectors 

IRFDF 

Subfile  device  number  for  barriers 

IT 

Alphanumeric  "T" 

ITY 

Vehicle  type 

J 

Index 
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Index 

LAST 

Alphanumeric  "L" 

NB 

Number  of  barriers 

NG 

Number  of  absorptive  ground  strips 

NQ 

Number  of  vehicle  types 

NQC 

Number  of  groups  within  one  vehicle_  type 

NQS 

Vector  notation  for  number  of  vehicle  groups 

NR 

Number  of  roadvjays 

NRC 

Number  of  receivers 

NSEC 

Section  number 

NSMl 

NSEC-1 

RDIN 

Vector  notation  for  initialization  parameters 

VALUE 

Initialization  parameter 

VEH 

Vehicles  per  hour 

VEXPH 

Vehicles  per  foot,  vehicles  per  hour 

X 

Coordinate 

XGl 

X-coordinate  for  end  point  of  ground  strip  center  line 

XlvIH 

Speed  in  mph  for  one  group  of  vehicles 

XMPH 

Vector  notation  for  vehicle  speed 

XNIGHT 

Height  adjustment  for  receiver  coordinate 

XRC 

X-coordlnate  of  receiver 

XT 

X-coordinate  for  end  point  of  road  section 

Y 

Coordinate 
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YGl 

YRC 

YT 

Z 

ZGl 

ZRC 

ZT 


Y-coordinate 

Y-coordinat e 

Y-coordinate 

Coordinate 

Z-coordinate 

Z-coordinate 

Z-coordinate 


for  end  point  of  ground  strip  center  line 
of  receiver 

for  end  point  of  road  section 


for  end  point  of  ground  strip  center  line 
of  receiver 

for  end  point  of  road  section 


3.4  SubroutineGEOMRY 
3.4.1  Overvi ew 

The  bulk  of  the  computational  effort  of  this  subroutine  is 
devoted  to  finding  segments  of  a straight  road  section  for 
which  sound  from  all  points  along  the  segment  reaching  a'  par- 
ticular receiver  is  nearly  equally  attenuated  by  barrier  dif- 
fraction or  ground  absorption  or  is  equally  amplified  by 
reflections . 


Because  of  the  complexity  of  this  subroutine,  a detailed  de- 
scription with  a series  of  flow  diagrams  is  presented  in 
Sec.  3-4.2.  Before  reviewing  these  diagrams,  however,  the 
reader  should  become  familiarized  with  the  basic  philosophy 
of  the  computational  scheme  outlined  below. 

As  shown  in  the  following  diagram  giving  an  overview  of  sub- 
routine GEOMRY,  checks  are  performed  on  two  different  levels  — 
first  on  all  the  barriers,  then  on  all  the  absorptive  ground 
strips,  and  finally  on  all  the  reflectors.  On  the  lower  level 
preliminary  checks  are  made  mostly  in  the  two  dimensions  of 
the  x-y  plane  to  determine  whether  a barrier,  a ground  strip. 
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AN  OVERVIEW  OF  SUBROUTINE  GEOMRY 


42 


1 


I 

I 

m 

1 

I 

1 


1 

1 

i 

1 

1 

L 


u 

f 


or  a reflector  can  be  at  all  effective  for  the  road  segment 
and  the  receiver  under  consideration.  During  these  checks, 
subfiles  are  created  for  a smaller  number  of  barriers,  ground 
strips,  and  reflectors  that  need  to  be  Investigated  more 
closely  at  the  higher  level. 


At  this  higher  level,  break  points  between  road  segments  pos- 
sibly affected  and  unaffected  by  barriers,  ground  strips,  or 
reflectors  are  determined  from  checks  first'  made  in  the  two 
dimensions  of  the  x-y  plane.  If  the  road  segment  between 
the  initial  point  and  the  break  point  is  unaffected,  then 
no  further  checks  in  three  dimensions  are  made  on  the  actual 
effect  of  the  remaining  road  segment,  but  the  break  point  is 
considered  as  the  new  end  point  of  the  road  segment.  The 
remaining  road  segment  is  investigated  after  the  calculations 
for  the  first  road  segment  are  finished. 


If  the  road  segment  between  the  initial  point  and  the  break 
point  is  possibly  affected  by  a barrier,  a ground  strip,  or 
a reflector,  the  break  point  is  taken  as  a preliminary  end 
point  until  the  subsequent  analysis  in  three  dimension  shows 
whether  there  is  an  actual  effect  that  makes  it  necessary  to 
keep  the  break  point  as  an  end  point  or  whether  there  is  no 
effect  and  the  break  point  can  be. dropped. 


Since  all  the  checks  are  performed  first  for  diffraction,  then 
for  absorptive  ground  strips,  and  finally  for  reflections,  the 
Initial  road  section  may  be  reduced  several  times  by  each  group 
of  elements.  The  procedure  of  subdividing  a road  section  is 
diagrammed  in  Fig.  10.  The  final  reduction  in  each  group  is 
stored  in  terms  of  end  points  XR2D  for  barriers,  XR2G  for 
absorptive  ground  strips,  and  XR2  for  reflectors.  By  this 
means,  duplicate  calculations  are  avoided.  After  the  noise 
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CONSTANT  ATTENUATION  BY  DIFFRACTION 
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UNIFORM  CONTRIBUTION  OF  REFLECTED  RAYS 
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FIG. 10  SUBDIVISION  OF  ROAD  SECTION 
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from  the  road  segment  XR10,  XR2  is  analyzed,  the  computation 
continues  for  possible  reflections  from  the  road  segment  XR2, 
XR2G.  'After  the  noise  from  the  road  segment  XRl , XR2G  is 
completely  analyzed,  the  bomputation  continues  with  considera- 
tion of  ground  absorption  and  reflections  for  the  road  segment 
XR2G,  XR2D.  The  designation  of  end  points  during  the  course 
of  this  analysis  is  illustrated  in  Fig.  11.  Finally,  the 
computation  continues  with  the  entire  problem  of  barrier  dif- 
fraction, ground  absorption,  and  reflection  for  the  remaining 
road  segment  XR2D,  XR20  until  the  noise  from  the  entire  road 
section  XR10,  XR20  has  been  analyzed.  At  that  point,  program 
control  is  returned  to  the  main  program. 

3.4.2  Detailed  flow  diagrams 

The  flow  diagrams  presented  at  the  end  of  this  section  show  the 
operations  and  important  branch  points  of  the  computations  that 
are  indicated  in  the  overview  diagram  presented  above.  Again, 
statement  numbers  are  shown  so  that  reference  may  be  made  to 
the  corresponding  block  of  code  in  the  program  listing  given 
in  Sec.  3.6.  The  following  comments  are  provided  to  supplement 
the  information  contained  in  the  flow  diagrams. 

» 

PTetimi-naTy  oomputati-ons 

The  main  computation  shown  in  the  first  flow  chart  concerns  the 
distance  between  the  near  point  XNPT  on  the  source  line  and  the 
receiver  XRC.  This  distance  enters  into  the  formula  for  the 
mean  energy  level  and  is  a quantity  that  is  characteristic  of 
the  entire  road  section. 
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FIG. 11  FURTHER  SUBDIVISION  OF  ROAD  SECTION 
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PreliminaTy  selection  of  diffracting  barriers 


The  computations  are  made  in  the  two  dimensions  of  the  x-y 
plane,  and  no  checks  are  made  for  the  height  of  barriers.  A 
barrier  whose  top  line  crosses  the  roadway  in  a projection  on 
the  x-y  plane  is  considered  illegal  and  causes  the  program  to 
abort.  Thus,  even  a barrier  crossing  underneath  an  elevated 
roadway  is  illegal.  A barrier  stored  on  file  #IDF  is  not 
necessarily  high  enough  actually  to  diffract  sound  rays 
travelling  from  the  road  section  to  the  receiver.  The  storage 
of  the  relative  location  (see  SUBROUTINE  BLOCN)  saves  consider- 
able computation  for  very  high  barriers ' that  shield  the  entire 
road  section  against  the  receiver. 


Preliminary  selection  of  ah  sorbing  ground  strips 

The  computations  are  entirely  analogous  to  the  preliminary 
selection  of  diffracting  barriers.  No  check  is  made  for  the 
distance  of  the  ground  strip  center  line  from  the  road  section. 
Therefore,  a road  section,  or  a part  of  it,  which  is  inadver- 
tently specified  on  a ground  strip  will  not  be  detected  in  all 
cases  by  the  program. 

Diffraction  of  direct  ray 

The  initialization  of  a path  length  DELP0  = -0.2  means  that  the 
attenuation  by  diffraction  is  set  equal- to  zero  in  the  frequency 
band  centered  at  500  Hz. 

The  test  for  whether  or  not  a barrier  possibly  shields  the 
entire  road  is  made  in  the  x-y  plane  and  doe.s  not  account  for 
the  actual  height  of  the  barrier. 
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The  height  difference  for  sound  rays  from  cars  above  the  barrier 
provides  a rather  conservative  check  for  possible  effects  of 
diffraction.  A positive  height  difference  means  the  source 
is  visible  from  the  receiver  location. 

Only  the  most  effective  of  a number  of  parallel  barriers  is 
considered;  thus,  the  barrier  having  the  strongest  diffraction 
is  retained.  Note  that  this  barrier  is  not  necessarily  the 
highest  one.  It  is  just  as  likely  to  be  the  barrier  nearest 
to  the  road  or  to  the  receiver . 

The  calculation  of  a road  segment  with  sufficiently  similar 
diffraction  is  performed  for  cars  only.  It  is  then  assumed 
that  the  same  result  holds  for  trucks  having  a different 
source  height. 

The  check  for  maximum  diffraction  of  sound  in  the  frequency 
band  centered  at  500  Hz  is  made  for  trucks  only,  because,  the 
higher  source  height  of  trucks  will  result  in  smaller  dif- 
fraction than  for  cars.  - If  this  test  reveals  a path  length 
difference  of  more  than  12.5  ft,  further  checks  on  other 
barriers  can  be  skipped  for  the  road  segment  considered. 

If  a very  high  barrier  .shields  the  entire  remaining  road 
section,  no  further  tests  for  diffraction  of  the  direct  ray 
are  made.  This  procedure  may  result  in  some  inaccuracies 
for  very  long  road  sections  for  which  the  diffraction  of 
sound  from  points  near  to  the  receiver  is  very  strong  but 
for  which  the  diffraction  of  sound  from  remote  points  might 
be  considerably  weaker.  However,  these  Inaccuracies  are  no 
more  serious  than  those  resulting  from  the  assumption  of  the 
arbitrary  value  of  12.5  ft  for  the  maximum  path  length  dif- 
ference, and  are  thus  neglected. 
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Ground  abaorption 

Absorptive  ground  strips  are  disregarded  for  sound  from  a road 
segment  whenever  an  attenuation  different  from  zero  is  found 
due  to  barrier  diffraction.  The  criterion  for  this  decision 
may  be  modified  in  the  future  to  be  5 or  10  dB , corresponding 
to  path  length  differences  of  0 or  0.1  ft  at  500  Hz,  if  field 
experience  shows  that  this  is  necessary.  Furthermore,  the 
decision  to  neglect  ground  strips,  which  is’  related  to  certain 
heights  of  rays  above  the  center  line  of  the  ground  strip, 
is  not  based  on  experimental  evidence  and  may  require  modi- 
fications. 

Preliminary  selection  of  reflectors 

By  the  introduction  of  an  image  receiver  location,  the  reflec- 
tion problem  becomes  similar  to  the  diffraction  problem.  A 
reflector  in  the  path  from  the  road  segment  to  an  image  receiver 
is  effective  whenever  a barrier  in  the  path  from  the  road  seg- 
ment to  a receiver  strongly  diffracts  the  sound. 

In  addition  to  the  preliminary  checks  made  above  for  diffracting 
barriers,  the  intensity  of  the  direct  (or  diffracted)  sound 
from  the  road  segment  considered  is  compared  to  the  potential 
maximum  contribution  of  each  reflection.  This  check  assures 
that  only  essential  reflectors  are  considered.  Since  a re- 
flector might  not  be  high  enough  or  reflections  might  be 
.strongly  attenuated  by  diffraction  at  additional  barriers 
in  the  ray  path,  the  reflectors  found  at  this  stage  are  con- 
sidered "potential"  reflectors. 
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Also  involved  in  the  part  entitled  "Preliminary  Selection  of 
Reflectors"  is  the  preliminary  selection  of  barriers  which 
can  possibly  diffract  the  reflected  sound.  ’ The  check  per- 
formed is  based  on  the  length  and  orientation  of  a barrier 
section  and  on  the  distance  from  the  image  receiver  relative 
to  the  respective  parameters  of  the  road  section  and  the 
receiver  considered. 

Catcutation  of  potential  reflections 

The  computations  in  this  section  are  analogous  to  those  made 
for  the  diffraction  problem,  except  for  a simplified  decision 
concerning  the  height  of  the  reflector.  If  the  height  dif- 
ference above  the  reflector  of  a sound  ray  from  the  nearest 
point  on  the  road  segment  leads  to  the  acceptance  or  rejection 
of  the  reflector,  no  further  check  is  made  for  other  source 
points  on  the  road  segment. 

Since  diffraction  of  the  reflected  ray  is  not  yet  checked,  a 
reflection  is  still  called  a potential  one  and  the  end  point 
of  the  road  segment  is  preliminary. 

Consideration  of  diffraction  before  reflection 

Checks  for  barriers  in  the  area  defined  by  projections  of  the  • 
four  points  XRl,  XR2,  XRB3,  and  XRB4  into  the  x-y  plane  are 
made  by  considering  first  the  triangle  containing  image  rays 
and  then  the  triangle  formed  by  the  road  segment  and  the  image 
receiver . 
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If  a barrier  is  found  to  be  high  enough  for  possible  diffraction, 
the  reflector  is  checked  again  to  determine  whether  or  not  it 
is  high  enough  to  reflect  the  diffracted  rays  (which  now  come 
from  an  effective  source  that  might  be  considerably  higher 
than  the  roadway). 

Calculation  of  the  path  length  difference  between  diffracted 
and  direct  rays  from  only  the  near  point  of  the  road  segment 
implies  that  the  diffraction  of  sound  rays  from  other  source 
points  is  about  the  same. 

Very  strongly  diffracted  reflections  are  neglegted.  The  deci- 
sion is  made  on  the  basis  of  the  diffraction  of  sound  from 
trucks,  since  rays  from  cars  are  even  more  strongly  diffracted. 

Consideration  of  diffraction  after  reflection 

Checks  are  made  for  diffracting  barriers  in  the  triangle. defined 
by  projections  of  the  reflector  segment  XRB3,  XRB4  and  the  re- 
ceiver XRC  onto  the  x-y- plane. 

After  a barrier  has  been  found  which  is  high  enough  for  possible 
diffraction,  the  reflector  is  checked  to  determine  whether  or 
not  it  is  high  enough  to  reflect  sound  towards  the  top  line 
of  the  diffracting  barrier,  which  might  be  considerably  higher 
than  the  receiver. 

The  calculation  of  the  path  length  difference  implies  simplify- 
ing assumptions  similar  to  those  for  the  problem  of  diffraction 
before  reflection. 

Reflections  are  neglected  in  the  case  of  diffraction  before  and 
after  reflection  and  in  the  case  of  very  strong  diffraction 
after  reflection. 


51 


Sit 


f." 


y 


».in  ,l:t 


■r-,.;  'vv>  ivH’ T0^^«ra4p’^‘ 

* i i ij'  -.I  /y?  ftj  •■/  fv^' 


,w^ 


I 


FvnaZ  computat-Cons  and  summation  of  sound  intensities 

In  order  to  facilitate  future  calculation  of  higher  order  sta- 
tistical parameters  of  the  noise,  the  contributions  from  all 
reflections  are  used  to  calculate  a gain  factor  of  the  direct 
sound  from  a road  segment.  The  subsequent  computations  are 
directly  oriented  at  the  formula  of  Eq . (8)  for  the  mean 

sound  intensity. 

Also  intended  for  future  use  is  the  separate  computation  of 
angles  ANGl  and  ANG2,  of  which  only  the  absolute  value  of  the 
difference  is  employed  in  Eq.  (8). 

After  the  calculation  for  a road  segment  has  been  completed, 
a point  1 ft  beyond  the  end  point  of  the  segment  is  taken  as 
the  initial  point  of  the  next  road  segment.  This  procedure 
avoids  computational  complications  without  unduly  sacrificing 
accuracy . 
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Preliminary  Computations 


Find 

the 

point  XNPT  on  source  line  containing 
road  section  nearest  to  receiver  XRC 

Compute  angle  ANGl  subtended  at  receiver 
XRC  by  line  segment  XRl,  XNPT 
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Preliminary  Selection  of  Diffracting  Barriers 


For  all  sections  of  all  barriers  in  file 


Preliminary  Selection  of  Absorbing  Ground  Strips 


For  all  ground  strips  on  file  #IGA: 


Does  center  line  of  the  ground 
strip  cross  the  road  section? 


RETURN  ^ 


N 


Does  the  center  line  of  the 
ground  strip  or  a .part  of  it 
lie  in  the  triangle  formed  by 
the  road  section  XR10,  XR20 
and  the  receiver  XRC? 


1 r w 


Store  the  information  about 
the  ground  strip  and  the 
relative  location  on 
file  #IGTA 


Disregard  the 


ground 

strio 
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Diffraction  of  Direct  Ray 
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Ground  Absorption 
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Preliminary  Selection  of  Reflectors 


For  all  sections  of  all  barriers  on  file  #IBR: 
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Calculation  of  Potential  Reflections 
For  all  reflectors  on  file  i^^lREF: 


6o 


c: 

r: 


c 

( 

1 

1 

1 


Consideration  of  Diffraction  Before  Reflection 
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Consideration  of  Diffraction  After  Reflection 
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Final  Computations  and  Summation  of  Sound  Intensities 


65 


[ 

[ 

[ 

[ 

[ 

[ 

[ 

[ 

[ 

r 

r 

r 

r 

r 

r 

r 

r 

r 


66 


3.4.3  List  of  variables  for  subroutine  GEOMRY 


A 

ADST 
ANG 
ANGl 
ANG  2 
ANGI 
ANG IMG 
CQ 

DELM 

DEEP 

DELP0 

DELPl 

DEEP  2 

DEEPA 

DEER 

DEETA 

DIST 


Atmospheric  attenuation  factor,  ground  absorption 
parameter 

ANG/DIST 

Angle  subtended  at  receiver  by  road  segment 
Angle  in  Fig.  1 

Angle  in  Fig.  1 

Angle  subtended  at  image  receiver  by  road  segment 

Angle  subtended  at  image  receiver  by  barrier  section 

Factor  accounting  for  standard  deviation  of  reference 
level 

Mean  path  length  difference 
Path  length  difference 

Maximum  path  length  difference  for  diffraction  of 
direct  ray 

Maximum  path  length  difference  for  diffraction  before 
ref le  ction 

Maximum  path  length  difference  for  diffraction  after 
reflection 

Path  length  difference  for  diffraction  of  direct  ray 
Path  length  difference  for  reflected  ray 
Distance  along  the  roadway 

Distance  from  the  receiver  to  the  source  line 
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DISTI 

Distance 

from  the 

image 

receiver 

to 

the  source 

line 

DISTJ 

Distance 

from  the 

image 

receiver 

to 

the  diffracting 

barrier 

DL 

Mean  path 

, length 

over  an 

1 absorptive 

ground  strip 

DNl 

Distance 

from  the 

receiver  to  the 

ne 

arest  point 

of 

the  road 

segment 

DNII 

Distance 

from  the 

image 

receiver 

to 

the  nearest 

point 

of  the  road  segment 

DN2 

Distance 

from  the 

image 

re  ceiver 

to 

the  nearest 

point 

of  the  diffracting  barrier 

DRl  Distance  from  the  receiver  to  the  initial  point  of 


the  road  segment 

DRK  Distance  from  the  receiver  to  the  preliminary  end 
point  of  the  road  segment 

FB  Attenuation  factor  accounting  for  diffraction  and 
reflections 

FCTR  Weighting  factor  for  reflections 
FG  Ground  attenuation  factor 

HDIFA  Height  of  ray  from  source  point  XRl  to  receiver  XRC 
above  barrier 

HDIFF  Height  of  ray  above  barrier,  reflector,  or  ground 
strip 

HGA  Data  for  effective  height  of  ground  cover 
I Index 
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lA 

Alphanumeric  "A" 

I BAR 

Barrier  number 

IBR 

File  device  number  for  barriers 

ICODE 

Number  for  intermediate  print  out 

IDF 

Subfile  device  number  for  barriers 

IDXR 

Number  of  reflections 

lERR 

Error  index 

IGA 

File  device  number  for  absorptive  ground  strips 

IGRA 

Ground  strip  number 

IGTA 

Subfile  name  for  absorptive  ground  strips 

IF 

Index  for  frequency  bands 

II 

Index 

IK 

Frequency  band  number 

IKIND 

Index  for  kind  of  absorptive  ground  cover 

IP 

Frequency  band  number 

IQ 

Index 

IR 

Alphanumeric  "R" 

IRD 

File  device  number  for  roadways 

IREF 

File  device  number  for  reflectors 

IRFDF 

Subfile  device  number  for  barriers 

ISEG 

Barrier  section  number 

c 


ITRIG 


KAR 

KBAR 

KBARl 

KBAR2 

KCD 

KDIFF 

KF 

KGA 

KIMG 

KREF 

KRFDF 

KTRIG 

LOG 

MDIFF 

MODD 

NB 

NBSEC 

NDIFF 

NF 

NG 

NIMG 


Trigger 

Alphanumeric  indicator  for  type  of  barrier 

Total  number  of  barrier  sections;  index  for  barriers 

Reflector  number  on  subfile  IREF 

Diffractor  number  on  subfile  IRFDF 

Indicator  for  relative  location  of  b-arrier 

Barrier  number  on  subfile  IDF 

Index  for  frequency  bands 

Number  of  ground  strips  on  file  IGTA 

Reflection  number 

Reflector  number  on  subfile  IREF 

Barrier  number  on  subfile  IRFDF 

Indicator  for  intersection  of  barrier  or  ground  strip 
Indicator  for  relative  location  of  barrier_^or  ground 
strip 

Indicator  for  diffraction  before  reflection 

Indicator  for  diffraction  of  direct  ray 

Number  of  barriers 

Number  of  barrier  sections 

Number  of  barriers  on  subfile  IDF 

Number  of  frequency  bands 

Number  of  absorptive  ground  strips 

Number  of  reflections 
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NQ 

NQQ 

NQS 

NR 

NREF 

NRFDF 

PP 

RATIO 

RDIN 

T1 

T2 

VEXPH 

XBl 

XB2 

XDBl 

XDB2 

XDB3 

XDB^ 

XGl 

XG2 

XG3 

XG^ 


Number  of  vehicle  types 

Number  of  groups  within  one  vehicle  type 
Vector  notation  for  number  of  vehicle  groups 
Number  of  roadways 

Number  of  reflectors  on  subfile  IREF 
Number  of  barriers  on  subfile  IRFDF  . 

Frequency  band  number 

V/eighting  factor  for  reflected  rays 

Vector  notation  for  initialization  parameters 

Temporary  variable 

Temporary  variable 

Vehicles  per  foot 

Initial  point  of  barrier  on  subfile  IDF 

End  point  of  barrier  on  subfile  IDF 

Initial  point  of  barrier  on  subfile  IRFDF 

End  point  of  barrier  on  subfile  IRFDF 

Initial  point  of  effective  barrier  segment 

End  point  of  effective  barrier  segment 

Initial  point  of  center  line  of  absorptive  ground 

strip 

End  point  of  center  line  of  absorptive  ground'  strip 
Initial  point  of  effective  ground  strip  segment 
End  point  of  effective  ground  strip  segment 
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XIMG  Vector  of  image  receivers  for  all  reflections 

XJ  Preliminary  end  point  of  effective  reflector  segment 
XK  Preliminary  end  point  of  road  segment 
XKI  Preliminary  end  point  of  image  road  segment 
XLE  Mean  intensity 

XLREF  Vector  notation  for  reference  intensities 

XNl  Point  on  road  segment  nearest  to  receiver 

XNII  Point  on  road  segment  nearest  to  image  receiver; 
point  on  image  road  segment  nearest  receiver 

XN2  Point  on  barrier  segment  nearest  to  image  receiver 

XNPT  Point  on  source  line  nearest  to  receiver 

XNPTI  Point  on  source  line  nearest  to  image  receiver; 

Oil  jLiua.^0  oOvaiTcc  Tunc  dgOjCcco  ncccivcn 

XR  X-coordinate  of  receiver 
XRl  Initial  point  of  road  segment 
XRL0  Initial  point  of  road  section 
XRII  Initial  point  of  image  road  segment 
XR2  End  point  of  road  segment 
XR20  End  point  of  road  section 

XR2D  End  point  of  road  segment  vjith  constant  attenuation 
by  diffraction 

XR2G  End  point  of  road  segment  with  constant  attenuation 
by  ground  absorption 

XR2I  End  point  of  image  road  segment 
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XRBl  Initial  point  of  reflector  on  subfile  IREP 
XRB2  End  point  of  reflector  on  subfile  IREP 
XRB3  Initial  point  of  effective  reflector  segment 
XRB4  F’''d  point  of  effective  reflector  segment 
XRC  Receiver  point 
XRCI  Image  receiver  point 
YR  Y-coordlnate  of  receiver 
ZN10  Z-coordinate  of  XNl  or  XNII 
ZR  Z-coordinate  of  receiver 
ZS  Height  adjustmient  for  vehicles 

3.5  Description  of  Other  Subroutines  and  Functions 
SUBROUTINE  FILES 

FILES  defines  file  device  numbers  for  the  three  files  with 
roadway,  barrier,  and  absorptive  ground  strip  data" created 
by  the  subroutine  INPUT  and  for  the  four  subfiles  with  barrier 
and  ground  strip  data  created  by  the  subroutine  GEOMRY. 

SUBROUTINE  INTER  (IQ) 

INTER  interpolates  between  stored  A-weighted  overall  and  octave 
band  reference  levels  to  find  reference  intensities  in  fi-e- 
quency  band  numbers  IQ  as  a function  of  the  average  speed  of 
vehicles  in  a group. 

FUNCTION  AMAG  (XIV,  X2V) 

AJ-IAG  is  a function  subroutine  v/hlch  calculates  the  distance 
betv.'een  the  points  XIV  and  X2V. 
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FUNCTION  ANGLE  (XIV,  X2V,  X3V) 

ANGLE  is  a function  subroutine  which  calculates  the  angle  at 
point  X3V  that  is  subtended  by  the  line  between  XIV  and  X2V. 

The  angle  is  positive  with  values  between  0 and  tt  . 

FUNCTION  BARFAC  (KF,  DEEP) 

BARFAC  is  a function  subroutine  which  calculates  the  attenuation 
factor  due  to  barrier  diffraction  as  a function  of  the  frequency 
band  number  KF  and  of  the  path  length  difference  DEEP  of  the 
diffracted  ray  path  relative  to  the  direct  ray  path. 

SUBROUTINE  BLOCN  (XIV,  X2V,  X3V,  X4V,  X5V,  X6V,  LOC) 

BLOCN  calculates  the  location  LOC  of  a barrier  or  a ground 
strip,  defined  by  the  end  points  X4V  and  X5V,  relative  to  the 
triangle  formed  by  two  points  XIV  and  X2V  on  a road  segment 
and  the  x-eceiver  point  XjV  (see  Fig.  12).  The  point  X6V  is 
identical  with  X2V  if  both  points  X4V  and  X5V  are  outside 
the  triangle.  If  one  of  the  points  X4V  or  X5V  is  inside  the 
triangle,  X6V  is  the  intercept  of  a line  through  X3V  and  X4V 
v;ith  the  road  segment  XIV,  X2V.  If  both  points  X4V  and  X5V 
are  inside  the  triangle,  then  the  point  X6V  is  the  intercept 
of  the  road  segment  XIV,  X2V  with  the  line  through  X3V  and 
either  X4V  or  X5V  depending  on  v/hich  intercept  is  nearer  to 
the  poirif-.  XIV. 

LOC  = 0 denotes  that  the  line  X4V,  X5V  is  outside  the  tri- 
angle XIV,  X2V,  X3V. 

LOC  = 1 denotes  that  the  line  X4V,  X5V  intersects  the  line 
XIV,  X3V. 

denotes  that  the  line  X4V,  X5V  intersects  the  line 
X2V,  X3V. 
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FIG. 12  NOTATION  FOR  SUBROUTINE  BLOCN 
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LOG  = 3 denotes  that  the  line  X4V,  X5V  intersects  both  lines 
XIV,  X3V  and  X2V,  X3V. 

LOG  = 4 denotes  that  the  line  X4V,  X5V  is  completely  inside 
the  triangle  XIV,  X2V,  X3V. 

FUNCTION  DEL  (XIV,  X2V,  X3V,  X4V,  HDIFF,  DNl) 

DEL  is  a function  subroutine  which  calculates  the  path  length 
difference  between  the  ray  emanating  from  the  source  point 
XIV  and  diffracted  over  the  barrier  X3V,  X4V  towards  the 
receiver  X2V  and  the  direct  ray  from  XIV  to  X2V.  HDIFF  is 
the  height  of  the  direct  ray  above  the  barrier,  and  DNl  is 
the  path  length  traveled  by  the  direct  ray. 

SUBROUTINE  ENDPT  (XIV,  X2V,  X3V,  a4V,  X5V,  X6V,  KTRIG,  lERR) 

ENDPT  calculates  the  nev;  end  point  of  a road  segment  from 
wnich  sound  is  unif ox'iiil^  attenuated  during  propagation  towards 
the  receiver  X3V.  The  original  road  segment  has  the  initial 
point  XIV  and  the  end  point  X2V.  A barrier,  reflector,  or 
ground  strip  has  the  end  points  X4V,  X5V.  X6V  is  the  break 
point  on  the  road  betv^een  segments  that  are  affected  and 
unaffected  by  the  barrier,  reflector,  or  ground  strip.  If 
the  entire  road  segment  XIV,  X2V  is  unaffected,  KTRIG  = 0 
and  X6V  is  identical  with  X2V.  If  the  entire  road  segment 
XIV,  X2’^^  -is  affected,  KTRIG  = 1 and  X6V  is  identical  with 
X2V.  If  a break  point  X6V  is  at  a distance  of  less  than 
0.51  ft  from  the  initial  point  XIV,  this  break  point  is  dis- 
regarde-^  If  the  initial  road  segment  XIV,  X6V  is  affected 
by  a barrier,  reflector,  or  ground  strip,  KTRIG  = 1,  the 
break  point  X6V  is  moved  tov/ards  XIV  by  0.5  ft,  and  the  new 
end  point  is  denoted  as  X6V.  If  the  initial  road  segment 
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is  unaffected,  KTRIG  = 0,  the  break  point  X6V  is  moved  towards 
XIV  by  0.5  ft,  and  the  new  end  point  is  denoted  as  X2V. 

The  error  indicator  lERR  is  set  in  the  SUBROUTINE  MOVE. 
FUNCTION  HEIGHT  (XIV,  X2V,  X3V,  X4V) 

HEIGHT  is  a function  subroutine  which  calculates  the  height 
of  a ray  from  point  XIV  to  point  X2V  above  the  line  with  the 
end  points  X3V,  X4V. 


FUNCTION  lEPS  (XIV,  X2V,  X3V,  X4V,  DELI) 


lEPS  is  a function  subroutine  which  calculates  whether  the 
path  length  difference  DELI  is  sufficiently  similar  to  the 
path  length  difference  between  the  ray  emanating  from  XIV 
and  diffracted  over  the  barrier  X3V,  X4V  towards  the  receiver 
X2V  and  the  direct  ray  from  XIV  to  X2V.  If  the  two  path 
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wise,  lEPS  = 0. 


SUBROUTINE  IMAGE  (XIV,  X2V,  X3V,  X4V) 

IMAGE  calculates  the  location  of  an  image  receiver  X4V  relative 
to  a receiver  X2V  and  a reflector  v;ith  the  end  points  XIV, 

X2V. 


SUBROUTINE  INTCPT  (XIV,  X2V,  X3V,  X4V , X5V) 

INTCPT  calculates  in  the  x-y  plane  the  intercept  X5V  of  two 
lines  specified  by  the  two  pairs  of  points  XIV,  X2V  and  X3V, 
X4V. 


77 


r: 

[. 

[- 

[■ 


■1 

K 


FUNCTION  KCUT  (XIV,  X2V,  X3V,  X4V) 


KCUT  Is  a function  subroutine  which  checks  whether  two  line  seg- 
ments with  the  end  points  XIV,  X2V  and  X3V,  X4V,  respectively, 
cross  in  the  x-y  plane.  If  the  two  line  segments  cross,  . 

KCUT  = 1-  otherwise,  KCUT  = 0. 

FUNCTION  KPOS  (XIV,  X2V,  X3V) 

KPOS  is  a function  subroutine  which  checks  whether  the  point  X3V 
lies  between  two  points  XIV  and  X2V  on  a line.  If  it  does, 

KPOS  = 1;  otherwise,  KPOS  = 0. 

SUBROUTINE  MIDP  (XIV,  X2V,  X3V) 

MIDP  calculates  the  center  point  X3V  between  two  points  XIV 
and  X2V. 

SUBROUTINE  NRPT  (XIV,  X2V,  X3V,  X4V,  DIST) 

NRPT  calculates  the  point  X4V,  on  a line  specified  by  the  points 
XIV,  X2V,  that  is  nearest  to  the  point  X3V.  It  als_o  calculates 
the  distance  DIST  betv/een  the  points  X3V  and  X4V.  If  the  dis- 
tance is  zero,  DIST  = 1. 

SUBROUTINE  NRl  (XIV,  X2V,  X3V,  X4V,  DIST,  X5V,  DNl) 

NRl  calculates  the  point  X5V,  on  a line  segment  between  the  end 
points  XIV,  X2V,  that  is  nearest  to  the  point  X3V.  It  also 
calculates  the  distance  DNl  between  the  points  X3V  and  X5V. 

The  subroutine  employs  as  input  parameters  the  distance  DIST 
between  the  receiver  and  the  nearest  point  X4V  on  the  entire 
line . 

SUBROUTINE  REPLCE  (XIV,  X2V) 

REPLCE  sets  the  vector  X2V  equal  to  the  vector  XIV. 
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SUBROUTINE  SECTN  (XIV,  X2V,  X3V,  X4V,  X5V,  X6V,  X7V) 

SECTN  calculates  the  segment  X6V,  XyV  of  the  line  between  X4V 
and  X5V  that  is  bounded  by  the  two  lines  from  X3V  to  XIV  and 
from  X3V  to  X2V. 

SUBROUTINE  TRI  (XIV,  X2V,  X3V,  X4V,  X5V,  KTRI ) 

TRI  determines  whether  or  not  the  point  X4V  lies  in  the  triangle 
formed  by  the  points  XIV,  X2V,  X3V.  If  X4V  is  in  the  triangle, 
KTRI  = 1;  otherwise,  KTRI  = 0.  TRI  also  calculates  the  inter- 
cept X5V  of  two  lines  specified  by  the  points  XIV,  X2V  and  X3V, 
X4V,  respectively. 

FUNCTION  ZCOR  (XIV,  X2V,  X3V) 

ZCOR  is  a function  subroutine  which  calculates  the  z-coordinate 
of  the  point  X3V  so  that  the  point' X3V  lies  on  the  line  speci- 
fied by  the  points  XIV  and  X2V. 

SUBROUTINE  MOVE  (XIV,  X2V,  X3V,  DELTA,  lERR) 

MOVE  calculates  the  point  X2V,  on  the  line  specified  by  the 
points  XIV  and  X2V,  that  is  at  a distance  of  DELTA  feet  from 
the  point  XIV  in  the  direction  away  from  X3V.  MOVE  gives  the 
error  message  lERR  = 4 if  the  two  points  XIV  and  X3V  are 
identical . 

FUNCTION  TAN  (X) 

TAN  is  a function  subroutine  which  calculates  the  tangent  of  X. 
FUNCTION  ACOS  (X) 

ACOS  is  a function  subroutine  which  calculates  the  arc  cos  of 
X betw'een  0 and  tt. 
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T 


J 

1 

JL 

I 

I 

I 

I 

i 

1 

1 

I 

1 

I 


3.6  Computer  Program  Listing 

TRAFFIC  NOlS£  PFFDICTION  MODEL 
DIMENSION  SIG(2) 

di::ension  xp  1 ( 3 ) , xr2  ( 3 ) 

COMMON /IN0U/I3 D, IBR, IDF, IREF, IFF DF, IGA , IGT A 
COM MON /BEK 2/ NQ 
C OMMC  N / I N p r 1 /F  D I N { 1 5 ) 

COMMON/IRPT2/NR, NB, NG 

CC;iMCK/IN?T3/XRC  ( 15)  ,IRC(  15)  ,ZRC(  15)  ,NRC 
COMMON/DPIV 1/XMPH ( 5 ) 

COi;MOfC/DPIV2/KC)S  ( 2 ) , NF 
C3MMON/D5IV 3/CQ ( 2 ) , VEXPH ( 5, 2 ) , XLE ( 9 ) 
ccmmon/gfo ;/ibae, iseg , igra 
INTEGLE  title  (3'‘) 

EOUIVALEN’CE  ( F.DIN  ( 3 ) , SIG  ( 1 ) ) 

F0  = 2 
NG  = 0 
Nb  = C 

CALL  SEARCH  { 1 , 6HINPUT  , 1 , :2  ) 

CALL  FILES 

6 CALL  S E A R C H ( L , 6 H S 0 A D S , I R D , 0 ) 

CALL  SEARCH ( a , 6H3APRIR, IBF  , 0 ) 

CALL  SEARCH ( 4 , 6KGHABS  ,lGAr0) 

C/..LL  SEARCH  ( 4 . 6HDEFLEC  , iDi  , 0 ) 

CnLL  SEA^CH(4, 6 H G RT A B S , IG T A , 0 ) 

CALL  SEARCH ( 4 , 6HREFLEC , IREF, 0 ) 

CALL  SEARCH ( 4 , 6 HRLFDFR , IRFDF,  0 ) 

PEAD(4,1005)  (LITL-.(I)  ,1=1,30) 

T ' ‘O  T rf.  r*  / 'If,.  ON 

ViRITE{  1,  1035)  (7ITL3(I)  , 1=1,30) 

CALL  INPUT 
NF=PDIN ( 2 ) 

HRJTE(  1,  10.2)  (TITLE (I) ,1=1,30) 

V:RITE  ( 1 , 10J3  ) 

DO  11  1=1,  :i  Q 

CQ(I)=EXP(  .5*  (SIG  (I)  *,23Z26)*  2) 

11  CuNTI.NUE 

DC  60  1=1, NRC 
DO  15  u = 1 , i ; F 
XLE ( J )=g. 

15  CONTINUE 

call  search  ( 1 , 6HR0=-.DS  ,IRD,0) 

DO  50  11=1, NR 

DO  20  10=1, NO 

CALL  BEAD ( IRD , NQQ ,1,0) 

CALL  ES:^n  ( IPD,  VEXPH  { , IQ  ) , 2 * NQO  , 0 ) 

CALL  READ ( IBDfXHPH ( 1 ) ,2*NQO,0) 

NQ3 ( 10 ) =MOQ 
IF{NOQ.EO.  )GG  TO  20 
CALL  INTER (IQ) 
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20 


r>D 


5C 


55 


COKTINUS 

CALL  READ(IRD,NES5C. ..0) 

CALL  REAL/viSD,XP.1,6,0) 

DO  40  i\’-1#MRSEC 

CALL  READ  (I?.D,XS2, 6.0) 

CALL  GF.OnRy  ( X DC  ( I ) , Y E C ( I ) . Z R C ( I ) . X E 1 . X R 2 . I L R R ) 

IF ( lEER , EQ. 1 ) GO  TO  65 

IF  { lERP.  . EC> . 2 ) GO  TO  7 0 

IF ( lERE . EQ . 3 ) GO  TO  66 

IF(IEKR.EQ.U)  GO  TO  60 

CALL  KEPLCE ( XR2, XR 1 ) 

C C N T I N U E 
COLTINUE 

CALL  SEARCH ( U, 6HB0ADS  .IRD.0) 

DO  55  0 = 1.  liF 

XLE  ( J ) =34  . + 1 0 , + ALOG  1 0 ( XL?:  ( 0 ) ) 

CGHTIi-iuE 

WRITE ( 1,  10J4)  I.XRC(I)  . YBC(I)  .ZEC(I) 
IF(NF.EQ.I)  GO  TO  56 
WRITE ( 1 . 200  1 ) 

V;RITF(1,20}2)  (XLE{II).II  = 2.NF) 

WRITF (1.2003) 

WRITE ( 1.  1007)  XLS(1) 

IF(NF.LE.6^  GO  TO  60 
X LA  = 0.0 
DO  5 9 0 = 2. N F 

XLA=XLA+10.T**(XLE(O)/10.) 

CGRT7  KM  IF 

XLA=10.0*ALOG  ^0(XLA) 

WRITE ( 1 . 2 00  4 ) 


W R I T F 

( 1 . 

1 0 0 7 ) 

6 0 

C 0 K T I 

NUE 

GO  TO 

5 

65 

W'  RITE 

( 1 > 

10 

^•■6) 

GO  TO 

5 

66 

W R I T F ( 1 . 

10 

o9) 

GO  TO 

5 

70 

WRITE 

( 1. 

10 

CO 

GO  TO 

5 

002 

FORMA 

T(//3 

A 2 

A 

i 

003 

FORMA 

T { 1 

<>  \ l 

RFC 

0 

04 

forma 

T(4Xp 

13 . 

0 

05 

FORM  A 

T ( 3 

2 ) 

✓ 

I 

30  6 

F OR  M A 

T ( 4 

3HILL 

N.IGRA 


I V 


0 6 


009 


20 


f/ 1 

L>  I 


GAL  BARRIER  INTERSECTS  RQADWAY  FOR  RE 
1 . 3h5£;^I2. 2X.  2HB#I2 . 2X.  3h3S//l2  ) 

F ORHAT ( 48HILL: GAL  GROUND  STRIP  INTERSECTS  ROADWAY  F 
1H?.i^I2. 3riRSAI2 . 2X,  4H£GS#I2  ) 

FORi-1  AT  ( 25HTOO  MANY  R EFL  E CTI 0 N S . HC  V . I 2 . 4 H K#.I2|.4H 
format ( /14X. 2H635X,  3H  2 5 a X . 3 H 2 50 4 X . 3 H D 0 0 3 X . 4 H 1 0 0 0 3 X 
X 4HU0003X. 4H800r ) 


C#12f  2X  p 2HR#I2 

OR  REC#I2.2X.2 

S # . 1 2 ) 

, 4 H 2 0 0 0 3 X . 
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'3 

i 

I 

■j 

■j 


2302  FORMAT  ( 1?>;,  8F7  . 1 ) 

2003  FORMAT ( 12X, 6HA (500 ) ) 

.007  FORMAT  ( 1 , F7  . 1 ) 

200U  ROFMAT ( 14X,  ^HA) 

2009  F0RMAT(///15X,24HTRAFFIC  NOISS  P R U I C TI  OiN // ) 
END 
i 1 

Subroutine  fi7.x^s 

C0;'I10N/IK0U/IRD,  IBR,  IDF,  IREF^IRFDF,  IGA,  IGTA 
IHD  = 2 
I:3E  = 3 
IDF-U 


$ 1 


lRFDF=b 
IGA  = 7 
igt:.  = 3 

RETURN 

END 

subroutine  input 

INTEGER  ALPHA (37) 

DlilENSION  NQS  ( 2 ) , VTXPH  ( 5 , 2 ) , X HP K ( o , 2 ) , 

1 X ( 11 ) , I ( 1 1) , E ( 11  ) 

CGIiilON/lNOU/IFD  , IBB  , IDF  , IREF  , IBFDF  , IGA,  IGTA 
CO.TMON /3LK2/KQ 
CO;IriO”/IN?T  1/FDIN  ( 15) 

co;;mon/in?t2/nr,ns,  NG 

C0::50N/IM?T3/X^C  ( 'i5),YRC(i5),ZPCf15),NRC 
EQUIV ALENCT  ( EDIN ( 1 ) , XNIGHT ) 

DATA 
DATA 

data 

DATA 
D A T A 


IA/2HA  / 

IG/2HG  / 

IR/2HR  / 

IT/2HT  / 

LA5T/2HL  / 

5 READ(4, 100^)IGO,I1,I2 

GO  TO  ( 1 0 , 2 . , 3--  , 36  , ^ 0 , 5 A'  ) , IGO 
C garbage  data. .. .program  IRlTlALIZrti 


10 
1 1 


’ION  pakameter: 

WitITE(  1,20.;2) 

^EAD{4,  1001)VALU£, TDK, HAST,  (ALP  HA(J.),I=1f30) 
DIN  ( IDN  ) -V  ALU7- 
KRITF  ( 1 , ze,' A 1 ) value,  IDN,  ( ALPHA  ( I)  , I-~  1 , 32  ) 


I F ( I L A S T . N E H A S T ) G 0 
GO  TO  5 


T 0 1 1 


C VEHICLE  DATA 
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20  NE=I1 

C.i'.LL  SEA:v>.;{  ( 2 , 6HR0ADS  .IRD,0) 
DO  28  J=1; NF 
NSEC=  1 

DO  21  K = 1 , 0 

21  N'QS(K)=0 

22  PtAD(4,1002)  V5H,XhK,ITi'i.ILAS? 
NQS  ( ITY  ) =NOS  ( I'i  Y ) + 1 

NQC^NQS { ITY  ) 

V£XPH  ( ITY  )=VEH/XMK/5260  , 
XYiPH(I3QC,ITY)=XMH 
IF(lLAST.NY.LAST)GO  TO  22 
DO  23  K=i,:iy 


NOC  = NQS  I K ) 

CALL  KRITE  ( I?.D  , KOC  ,1,3) 

CALL  WRITE  ( IKD,  VEXPH  ( ‘ , K ) , 2>^NQC  , 0 ) 
call  write  ( IRD,  XhrH  ( , K ) , 2 N QC  , 0 ) 

COWTlIiUE 
FOADV’AY  DATA  SECTIONS 

2w  read ( U , 1003 ) X ( NSEC ) , Y ( MSEC ) , Z ( MSEC ) , ILAST 


23 


25 


27 

28 


IF  (HAST.  SO.  LAST)  GO  TO  25 
NSEC  = MSEC+  •; 

GO  TO  2^4 

IF ( MSEC- 1 , ME. 0 ) GO  TO  26 
WFITK ( 1 , 20  i? ) 
call  exit 
i;  1 = r:  c - " 

CALL  WRITS'  IKD , NSX  1 , , 3 ) 

CALL  WRITS {IFD,X( 1 ) ,2,0)  . 

CALL  WRITE{IRD, Y ( 1 ) ,2,0) 

CALL  WRITE  ( IxRD  , Z ( 1 ) , 2 , 0 ) 

DO  2 7 I = 2 , N S E C 

CALL  WRITE ( I RD , X ( I ) , 2 , 0 ) 

CALL  WRITE ( IRD , Y ( I ) , 2 , j ) 

CALL  WRITE(IF.D,Z(I)  ,2,2) 

C C M T I N U E 
COMTIMUS 

CALL  S K A R C H ( 4 , 6 H R 0 A D S , I R D , 0 ) 
GO  TO  5 


barrier  data  sections 

30  N3=11 

IF ( NB . EQ , 0 ) GO  TO  5 

call  SEARCH ( 2, 6HBARRIH, IBR, 0 ) 

DC  3 5 J=1,13E 


3 1 


hSEC=  1 


read ( 4 , 1003 ) X ( MSEC ) 
IFdLAST.EQ.IA.OR.I 
NSEC-NSECh-  : 

GO  TO  31 


, Y ( MSEC ) , Z ( NSEC ) 
LA3T.S0.IR )GO  TO 


ILA  ST- 
32 
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I 


I 


I 


i 


I 


> 


'll 

V:  [■ 


' A 


r 


••s 


32  I?  ( NS-:c- •>.  NE  . 2)  GO  TC  33 
WRITE  ( 1 r 2fe;'  11  ) 

CALL  EXIi 

33  w's:; i = NSEC- 1 

CALL  WRITE  ( I3R,  NSi^l  , ,0) 

call  W rite  (I3E,  hast,  ' , ) 

CALL  ;iRlTE(IBr,X{  1)  ,2,0 
C ; L L W R IT  E ( I B E; , Y ( 1 ) , 2 o ) 
call  write { IBP, Z ( O ,2,0) 

DO  3ii  I = 2,1ISSC 

CALL  WRIT- ( IBP , X ( I ) ,2,0) 

CALL  WRITE (IBE , Y ( I ; , 2, 0 ) 

CALL  WRITE(IBI:,Z(I)  ,2,0) 

3L  CONTINUE 

35  CONTINUE 

CJ.LL  SEARCH  ( 4, 6H3ARRIR,IBR,0  ) 

GO  TO  5 

C ABSOLBIKG  GROUND  STRIPS 

36  N G ~ I 1 

IE ( KG . EQ . S ) GO  TO  5 

CALL  SEARCH { 2 , 6HGEA3S  ,IGA,0) 

DO  57  I=1,NG 

REaD(4O00^)XG1,YG1,ZG1,BG 
CALL  WRITE ( IGA , XG 1,2,0) 

C.- LL  write  ( IGA  , YG  1 , 2 , 0 ) 

CALL  WRITE ( IGA , ZG 1 , 2 , 0 ) 

CALL  WKITE(IGA,BG, 2,0) 
p r.  ;0  ' iJi , 1 0 ) X G 1 , Y G 1 . 7,  G 1 « I D I i M 

IF  ( IDUM  . EQ  . LG  ) IDUN:^  1 
IF  ( IDUF  . Ep  . IT  ) IDlIi'i:=2 
CALL  WRITE ( IGA, XG 1 ,2,0) 
call  VJRITr:  ( IGA  , Y G 1 , 2 »>  0 ) 

CALL  WRITE ( IGA , ZG 1 , 2 , 0 ) 

CALL  WRITE(IGA, IDUH, : , 0 ) 

37  CONTINUE 

CALL  SEARCH ( 4 , 6HGRA BS  ,IGA,0) 

G 0 TO  5 

C RECEIVED  DATA 
4 0 N R C = I 1 

DO  41  I-ONRC 

BEAD  ( U , 'I  Jjioj  ) XHC  ( I ) , Y RC  ( I ) , ZRC  ( I ) , IDUH 
0FC(I)=ZRC(I)  + 1''.0 
4 1 CO  K T I K iJ  £ 

GO  TC  5 

50  Call  search ( 1 , 6HR0ADS  ,IRD,0) 


DO 

65 

J = 

1 , 

NF. 

^ t 

v«  K 

ITE 

( 1 , 

20 

.>2)  J 

T>  r' 

55 

X = 

1 , 

MQ 

CA 

ll 

REA 

D{ 

lED  , 

ii  Q C , 

1 , 

O 

^ ii 

LL 

REA 

D( 

IRD, 

V " X ? 

H( 

, K 

) , 2 * N 0 C , 0 ) 

CA 

L L 

REA 

D( 

IP.D, 

XFiPH 

( 1 

,K  ) 

,2=^KQC,0) 

I F 

( KOC. E 

Q c 

.■  ) GO 

.TO 

55 

DO 

5 4 

X - 

O 

■AQC 

8H 


\_1 

\ i 

I I 


r 

1 


w 


I 


5 4 V E X F F ( I , K } = V E X P H ( I , K ) * X M P H { I , K ) 5 8 . 

WRITE  ( ^,2^^U)K,  (I,v:-XPHU,K)  ,XMPH(I,K)  ,I=1,NQC) 
55  CONTINUE 

CALL  READ  ( IF.D,  NSEC  , 1 , 0 ) 

CALL  READfIP.D,XT,2,:) 

CALL  READ  (IFD,  YT,  2,  .:  ) 

CALL  READ (IRD, ZT, 2 . E ) 

WHIT?:  (1,203) 

WRITE ( 1 , 20W5 ) XT , YT , ZT 
NSFC=NSEC+ 1 
DC>  60  I = 2,.iSSC 
CALL  READ (IHD. XT, 2, A ) 

CALL  R S A D ( I R D , Y T , 2 , C ) 

CALL  READ  ( IRD,  ZT,  2,  ::  ) 

WRITE ( 1 , 2026 ) I, XT, IT, ZT 
60  CONTINUE 
65  continue 

call  SEARCH  ( a , 6HEOADS  ,IP.D,0) 

IF  { NP  . EQ  . 0 ) C-0  TO  8. 

call  search ( 1 , 6HEAHEIR, IBR, 0 ) 

DO  75  J=1,NB 

CALL  READ  ( IBP.,  N5EC  , ^ , 3 ) 

CALL  READ(ISR,ILAST, • ^0) 

CALL  READ(IER,XT,2,T) 

CALL  READ  (lER,  YT,  2 , ■•  ) 

CALL  READ  ( IBP. , ZI  , 2 , ) 

WRITE  ( 1 , 20'^7  ) J , ILAST 
P I T ( 1 . ? ' 5 > X T‘ ; 'i  Z T 

NSEC-NSEC+ i 
DO  70  1=2, NSEC 
CALL  R E A D ( I B R , X T , 2 r 0 ) 
call  READ  ( IBR,  YT,  2,  >,.  ) 

CALL  READ (IBR, ZT, 2, 0) 

W R I TE  (1,2 0 J 6 ) I ..  XT  , Y T , Z T 
70  CONTINUE 
75  CONTINUE 

call  search (4, 6KB AHRIR, IBB, 0) 

0 IF(NG.EQ.0)  GO  TO  9? 

CALL  S E A R C H ( , 6 H G E A.  B 3 , X G A , 0 ) 

DO  8 5 1='*,  !G 

CALL  P.Eau  URA  r XG  1 , 2 , 2 ) 

CALL  READ ( IGA,YG1 , 2, 0 ) 

CALL  RPJAD  ( IGA,  ZG  1,2,0) 

CALL  REJiD  (IGA  , BG  , 2 , / ) 

CALL  RE. AD f IGA, XT, 2, 0 ) 

CALL  READ ( IGA, YT, 2, 0 ) 

CALL  READ (IGA, ZT, 2, 0 ) 

CALL  READ  ( IGA  , 1 , 0 ) 

IF ( IDUK . Eg.  1 ) IDUM  = IG 
IF ( IDUH . Sg. 2 ) IDUM=IT 

WRITE ( 1 , 2Z  , 2 ) I, IDUK, XG 1 , YG 1 , ZG 1 , BG , XT, YT, ZT 
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85 

CONTINUE 

CALL  SEARCH ( 4 , 5HGRABS 

,1GA,0) 

9 0 

WRITE  ( 1 , 8 ) 

DO  9 5 I = 1 , N P. C 

WRITE ( 1 .20  ;6)I,XRC(I)  , 

YRC  ( I ) , ZRC  ( i. 

9 5 

C C T I N U E 

i ETUEN 

10  00 

r OEM  AT ( 315 ) 

100^ 

? OFMAT(S10.0,I5,4X,A1, 

10X, 30A2 ) 

1002 

F ORMAT ( 2E 1 . . r , 15 , 5X , A 

) 

« 1 3 3 
13  3^ 
i 2333 
! 2331 


F0K!1A3(5S1,'.:,A1  ) 
FCr;MAT(iJE13.3) 
FORMAT  ( 33hPHOGKA;'i  I 
FOFMAF  ('IX,F12.5.I1F 


2 002 

FORMAT ( /9HROADW 

2 00  4 

FORMAT  ( 9HNUMFlP. 

1 , 1 5 X , 2 E 1 3 . 4 ) ) 

2005 

FCP.MAT  ( 5HNUHBTR 

2006 

* 0 R /i  A T ( 3 / 1 2 , 2 X 

2007 

format ( /9H3ARP1 

2 0 3 6 

format ( /BKRECEI 

15V  12X,  1tiZ) 

2010 

format ( 26HINSUF 

V '5 

O il.  I 


23 1 1 F CPMAT ( 29HINSUFFICI 
20  12  r 0 ?.  K A T { 1 7 H b S 0 R B I R G 
k^y  , 1 H Z 1 2X  , 5HKID?H/i^ 
FORilAF{22X,  18HS0URC 
r n 


2313 

ii 


S3BK0UTZNE  GFOnFY( 
DliiEi:  SIOX  XG  1 ( 3 ) . X 
Lir. ENSIGN 
FELi3 ( 2 ) . DFL? 1(2) 
XB1(3),XB2(3),XDB 
XRB1  (3)  .XR32(3)  ,-X 
XH  1(3) , XH2 {3) , XR 1 
XK ( 5 ) , XKT ( 3 ) r XJ ( 3 
X N 1 ( 3 ) , X N 2 ( 3 ) r X N 1 
Z S ( 2 ) 

CGiiNOli/INOU/IRD  , IB 
CORHOi: /h1k2/l'Q 
COl-;MON/IKPT1/RDIN  ( 


X F. , Y R , Z R , X R 1 0 . X T;  2 0 , 1 E H R ) 

G 2 ( 3 ) , X G 3 ( 3 ) . X G a { 3 ) 

,DELP2(2)  ,FB(9,.  2)  , DEEP.  (2,  10)  .FG(9)  .HGA(2)  , 
1(3) ,XDB2(3) ,XDL3{3) ,XDB4( 3)  , 

RB3  ( 3 ) . XRB4  { 3 ) ,■  XKC  ( 3 ) , XKCl  ( 3 ) , 

I ( 3 ) , XR2I  ( 3 ) , XP  10  ( 3 ) , XP.20  { 3 ) ^ 

) , XNPT  ( 3 ) , X N p:;  I ( 3 ) . X N PT  J ( 3 ) , 

I ( 3 ) , XliiG  ( 3 1 0 ) ..  XR  2 D ( 3 ) , X R 2G  ( 3 ) , 

F,IDF,IBEFfIRFDF.IGA,IGTA 


15) 

MMGN/INPT2/NR , NB, NG 


C0;;ii0’:/DRiV2/ii0S  ( 2 
COMMON/DPIV3/CQ ( 2 ) 
cc.;mon/gso  1/iba 
CCiil-iON/INTEF  1/XLRF 
equivalence  (RDIN( 
I A / 2 H A / 
IH/2KF  / 


} .NF 

, VEXPH(5,2),XLE(9) 
G , I G R A 
F ( 5, 9. 2 ) 

5)  ,ZS(  1 ) ) 


T q T- 


FA 


r.  j 


O 


NITIALIZATION  PARA0EIERS) 

, 5X, 30A2  ) 

^ j.  3 ) 

3X, 5HVEH/H6X , 3HKPH/5H  TYPE, 12, UH  VEH/(3X,I2 

H X ; 2 X , 1 H y 1 2 X , 1 h Z / 4 X , 1 H 1 , 2 X , 3 E 1 3 . 4 ) 

3.4) 

1 3 , 2 X , : H ( A 1 , 1 H ) 1 0 X , 1 9 H B A R R I E R COORD  IN  FT) 

4X,2  HRECEIVEP  COORD  IN  FT/7H  N U M B ER 5 X , 1 H X 1 2 X , 1 

ENT  ROAD  SECTIONS) 

ENT  BARPIEB  SECTIONS) 

STRIP  AI3,2X,  1H(A1,  1H)//5H  PT  #7 X , 1 HX  1 2X , 1 H Y 1 
X , i K 1 2 X , 4 E 1 3 . 4 / 4 X , 1 H 2 2 X , 3 E 1 3 . 4 ) 

E COORD  IN  FT) 


..J 


J 


3 


-J' 


} 

J 


86 


!r'  ■■'« 

: ^®l'i2' f -I I ‘^«i^  . 

" '"  ■.  <Sf.'Sj 

.<?■".  *4#  ?»:  i'-r 

X'M'  H ’f^'^  'f  Sfeto  L&i.d  m^nu  '^'^41 0.  i'ig)M|^ 

'V'V';-  ■!^J^'’'' ■'  'i  *'''  Y r ^'1  V'M'V'f  ‘‘  ■*-  ' >n 
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DAI  A KGA(  1)  ,HGA(2)/  10.^30,/ 

XKC{  1 )=XR 

x:-c(  2 )=y:. 

XFC  ( 3 ) =?.R 
1ERF,  = 0 

CALL  iJRPT  (XB:  1C,  XR2C,  XRC^XNPT,  DIST  ) 

A N G 1 - A N G L E ( X E 'I  0 # X N P T , X B C ) 

ICODE=  1 

V.  F.  I T E ( 1 , 1 0 j r ) I C 0 D E , X R ■ 0 , X H 2 0 , X H C , X N P T 
PRELIXIFARY  SELECTION  CF  BARPIERS 
N D I F F = 0 

IF(i\'B,EO.0)  GO  TO  16 

CALL  SEARCH  ( I , 6H3AB  RIK  , IBP. , 0 ) 

CALL  SEARCH {2, 6 H D FE LEC , I D F , 0 ) 

K B A r.  = 0 

DO  15  IBAR=1,IIB 

CALL  READ ( IbR, NBSECr  ^0) 

C.ALL  RSAD(TB?,KAP,  1,0) 

CALL  READ (I3R,XB 1,6,0) 

DC  15  ISEG= 1 , NE5EC 
CALL  READ (IEF,X32, 6,0) 

KBAR  = KBAR+  , 

ICODE=2 

WRITE (1,1 00 1 ) ICODE , KBAR , XB 1 , XB2 
IF(KCUT(XRiP,XP,20,Xb  , ,Xb2  ) .N£„  1 )GC  TO  10 
I E H E = 1 
RETURN 

1C1  ft  AL  HT.OCN  t XK  10  . XF.20  , XPC  , XD  1 , XB2,  XK,  LOG  ) 
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15  CONTINUE 

CALL  SEARCH (4, 6HBARRIR, 13B, 0 ) 
CALL  SEARCH ( A , 6HDEFLEC, IDE , 0 ) 
PEILlnlNARY  SELECTION  OF  STKIps 

16  KGA=0 

IF ( KG. EG. 0) GO  TO  20 
CALL  SEARCH ( 1 , 6HGRABS  ,IGA,0) 
CALL  SEARCH(2,6HGRTABS,IGTA,0) 
DO  18  IGRA=1,NG 
C 7.  L L R E A D ( 1 G A , X G 1 , 6 , 0 ) 

CALL  READ(IGA,BG,2, A ) 

CALL  READ ( IGA , XG2 , 6 , 0 ) 

CALL  READ(IGA,IKIND,  ,0) 
IC0DE=3 
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WRITE  ( 1 , 1 ) ICODS,  IGRA/  XG  1 , XG2 

IF  ( KCUT  ( XR  IP,  XP.20,  XG  , X'G  2 ) . NE  . 'i  ) GO  TO  17 

IE?.?.  = 2 

RETURN 

17  CALL  BIOCN{XS10,XK2C.,XRC,XG1,XG2,a^,LOC) 
IF ( LOG . EO . ^ ) GO  TO  18 
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DO  30 
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IT"IG=0 

LC  12  KDir E=  1 , NDIr F 

r.  H M.  R E A n ( ,I  D F . K 3 A R . 1 . 0 ) 

CALL  READ(TD?,KCD, 1r0) 

CALL  K E A D ( T Lf , XB 1 , 5 , 0 ) 

CALL  READ (IBF,XB2, 6,0) 

IC0DE=5 

V.’PITE  (1,100")  ICODE  , KBAR  , XR  1 , X R 2 , X B 1 , X B 2 
IF ( KCD . EQ , 3 ) GO  TO  « L 

CALL  r H D P I { X r ^ X F:  2 , X R c , X 3 1 , X B 2 , X K , K T R 1 G , I E R R ) 
IF (IERR*EQ. A ) RETURN 
IF (KTRIG.EO,0)GO  TO  73 
GO  TO  US 

U0  CALL  EEPLCE(XR2,XK) 


4 5 r)CDD  = 0 
ICODF=6 

WRITS ( 1,1031 )IC0DS, KCD, XR2,XK 
DO  63  10=1,30 

CALL  NR  1 { XR  1 , XK  , Xac,  XI'IPT,  DIST,  XW  1 , DN  1 ) 
Z A 1 0 = X ! • 1(3) 

XN1  (3)=ZN1  " + ZS(IO) 

HDIFF  = HElGliT  ( XN1  , XEC,  X5  1 , XB2  ) 

IF (IQ, NE.  1 ) GO  TO  53 
IF ( HDIFF  . GT . 2A  . ) GO  TO  7J 
50  DF  LP  = DEL ( XN 1 , XRC, XB 1 , X32 , HDIFF, DN  1 ) 
IF(D£LP.GT.LELPC(IQ) )G0  TO  53 
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52  f r.ODD  . EQ  . 1 ) GO  TO  65 
G : TO  ivi 

5 5 I ( 1 0 . N E . , ) GO  TO  59 
:>?,  1 = a’iag  ( xp.c,  X?  1 ) 

I-'  (A5S(0R  -D'G1).LI.  .)  GO  TO  54 

n .iFA  = HEIGHT  ( X E , X HC  , X b 1 , X B 2 ) 

DEJ.PA  = DEL  (XR1  ,X^C,XB1  ,XB2,HDIf  A,DR1  ) 

: T L :•!  - ( D E L ? A + D E L ? ) / 2 . 

I C 0 D r = 1 7 

VRTTE(  1 . lOTDICCDE.DBLP^DELPA 

" =■  ■ ( A B S ( D E L .?  A - D E L ? ) - ii  . 1 - D T M / 5 0 , * ( 1 . + D E L M ) ) . L T . 0 . } GO  TO  5 5 
CALL  0IDP  (XR%XN1,XK) 

ICCDF=7 

V ? I T E ( 1 , 1 0 . • T ) I C 0 C L . X R ■ ^ X li  1 , X K 
DEL?  = D S L P A 

54  X-  ( lEPS  ( XK  , X?C  . XB  1 , XB2  , DKT.p  ) . EQ  , 0 ) GO  TO  58 
CALL  .']ID?  (XR-,XK,XK) 

i^ooE=a 

R I r E ( 1 , 1 0 X ? ) T c CD  E , X R p X K 
GO  TO  54 

55  rs;<=  AM  AG  (XRCpXK) 

I?  ( ABS  ( DRK-DX  1 ) . LT  ♦ ) GO  TO  58 

56  I:-  (TEPS  (XK,X?CpXB1pXB2,DELP)  .EQ.0)  GO  TO  58 
CALL  MIDP ( XN 1 . XK , XX ) 


GO  TO  56 
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IF  ( DTLP . LE , DTLP0 ( 1 ) ) GO  TO  52 
5 9 D E L P 0 (IQ)  = L E L ? 

T C o D V = 1 0 

AAITE  ( 1 p 1 0'C7  ) ICODE  p DELPky  ( 1 ) , DELP0  ( 2 ) 
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COL TI HUE 

t^'DELP0(2) .LT. 12,5)GO  TO  65 
I i E I G - 1 

irfXC0.NS.3.0T.XCD.NE.2)G0  TO  65 
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C ICCI^£=13 

C ( 1 , 10  ) ICCDI , X?  ■ , XR2  . XH2D 

.;  V V DFLP.;  ( 1 ) .GI  . -0 . 2 ) GO  XO  78 
I;  { XGA  . EQ  . ) GO  TO  7 8 

C....LL  S-:AEC8(  1,6HGRTABS,IGT?v.0) 
r 'J  7 7 I G H A = 1 , K G A 
CALL  READ ( IGTA , LOG ,1,0) 

CALL  READ ( IGTA , XG 1 , 6 , 3 ) 
call  head  ( IGIA  , XG2, 6 , 2)  ) 

C.\i,  READ(xGrA,BG,2,3) 

C.,Li.  HEAD  ( IGTA  , IKIAD  , i , 3 ) 

C ICODE=14 

C '••.AITS  ( 1 , 1 2 : 1 ) ICODS  . LOG,  XG  i , XG2 

I’:'(LCC.S0.3)G0  TO  7 4 

r:  T T,  en'DPT  : XH1  Xil2,  XAC,  XG  1 , XG2,  XK,  KTRiG,  ISRR) 
I " f lEHR . £Q  . a ) HLTURN 
IE ( KTHIG , SO . C ) GO  TO  77 
GA  TO  75 

74  C.sLL  R3PLCE  ( XR2,  XK  ^ 

75  CALL  >iiri  (XT', XK,XRC,XRPT, LIST, XN1,DN1  ) 
HlIFF-HSIGHT ( XN " , XAC, XG 1 , XG2 ) 

If  (HDIFF.GT,KGA(IKIXD)  )GO  TO  77 

CALL  sect:;  (XR1  ,XX,XHC,XG  1,XG2,XG3,XG4  ) 

LL=1 .57/(1  ./BG+'-  ./AMAG  ( XG3,  XG4  ) ) 
f'u  7 6 IK=1,MF 
I r = I K 

IFflK,SQ. 1 )I?=5 


7 6 


r» 


77 


7 8 
C PR" 


IF  ( IKIND  . SO  . 1 ) A=  ( ■/- . 00  1b^?P-.0 ,0328  ) *DL 
IF  ( IKIND  . E-O  , 2 ) A = 2 . / ■»=  * ( PP/'3 . 0 ) / 13  10  . ^DL 
IF(A .GT.3,)  A =3,0 
f L.  (IP). tFG  (IP) /1C.*.i.A 
IF(FG(IP)  .LT, 1 ,E-3)FG(IP)  = 1 .E-3 
CALL  RSPLCI ( XK, XR2 ) 

COATINUE 
IC0DE=15 

WRITE (1,1 032 ) ICODE , FG ( 

C'OKTINUS 

call  search ( 4 , 6HGRTAES 
CALL  3EPLC£(XS2,XR2G) 
li'.IUARY  SSLECIIOIi  OF  REFLiXTORS 
FRSF  = 0 

IF(MB,EQ.0)GO  TO  91 
XR?DF=2 
KBAR=0  . 

CALL  SSA-RCrK  1 , 6 H 5 A B P.  I R , IE  R , 0 ) 

ChLL  search ( 2, 6H&EFLEC,IRLF,0) 
call  search ( 2 , 6HEEF0FR, XRFDF, 0 ) 

LO  85  IBAH=1,NB 
call  READdBF  . NBSEC,  ,0) 

CnLL  READ(IS?.,KAR,  1 ,0) 

CALL  READ ( IBP, XRB 1,6,0) 

35  I5EG=  1 , N3SEC 
call  READ{I&R,XRB2.,6,3) 

KD.AR  = KBAR-<-  i 


1) , FG ( 9 ) 
,IGIA,0) 
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c 

I C R.‘  D 

E=16 
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V R i 

F ( 1 , 1 0 " 1 ) 

ICODE, 

K 5 

AR 
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B ; , X R 3 
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c L L 

NP.PT  ( XP  : 
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I C i‘  D 

F=  1 7 

c 

*?  r 

F(  1,  '!0.v) 

I C 0 D E . 

XE 
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A G X 

M G = G L E ( 

X R 3 1 , X 
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FCTE 

= ( n I S T * A N 

GIFiG  ) /DI 

S T 

IF  ( r 

Cl?  . LT , I . 

1 ) GO  T 

0 

80 

K r V D 

? = nrfdf^-  ; 

C A L 5 

WHITE ( IRFDF, KE 

AE 

ti 

f 1 
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WRITE ( IT 

FDF, XR 

E . 

. & 

C A L L 

S I T E ( I R 
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.6 
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E=:  1 3 

S ( 1 , 1 0 1 ) 

ICODE , 

N R F D 

li-  ( K 

AR  . RE  , IP ) GO  TO 

8 0 

CnLL 

ELOCN ( XR 

'1  , X R 2 

X R 

Cl 

I?(  L 

OC  . SQ  , j ) G 

0 TO  6 

'6 

NRZF 

— R E F + '’ 

C A L L 

WRITS ( IR 

E F , K 3 A 

R, 

' t 

CALL 

WRITS { IR 

E F , X F.  B 

1 . 

6 , 

CALL 

WRITE ( IR 

.ir  , XPP 

2, 

6 , 

ICOD 

?=19 

c 

WRIT 

2(1. 10 J1 ) 

I C 0 D .7  , 

LR 

r~\ 
■oJ  C 

3 0 

C rt  L L 

HEPLCE (X 

RB2,  X*". 

3 

) 

35 

cgn  r 

I N U E 

CALL 

SEARCH ( a 

, 6 HREF 

DF 

R, 

C 1 

^SG 

I « t: 

h'  L E •-  U r 

. O jL  1 1 

CALL 

3 S A ?.  C ri  ( 4 

, 6HBh? 

PI 

‘ 0 

C..  LL 

3 E A R C i1  ( 4 

, 6 H R F F 

L-  L 

f 

Q 0 

r T 

W /V  ^ AJ 

S E A R C H ( 1 

, 6HP.EF 

.W  Ij 

G 1 

I D X ? 

— .7 

IM  N 

RFF.EQ.OGO  TC 

179 

DO  175  KREF=1 

, NREF 

C /-» i_j 

READ ( IPE 

?, X3AR 

1 f 

• / 

CALL 

READ ( TPS 

? , X R B 1 

.6 

P *J 

C.;LL 

R S A D ( I R L 

F , X P 3 2 

. 6 

•'A 

? o 

c 

ICOD 

F;:^20 

n 

w PIT 

E ( 1 , 1 00  1 ) 

ICODF, 

KB 

AR 

CALL 

I;lAGE{  XR 

S 1 , X P.  B 

a , 

XP, 

CALL 

ENDPT ( XR 

1,XF2, 

XR 

Cl 

IF(I 

erh.?:q.4) 

RET  UR 

u 

lC0DE--2'i 


C WRITE ( 1 , 1 0j 1 ) ICODE , KTKI 

IF(KTRIG.EO,0)GO  TO  75 
CkLL  N'EPT  ( XR  1 ^ XE2  # XF.CI  , 
call  nr  1 ( XB 1 , XK , XRCI , XN 
C ICODE=22 


X ? 1 , X R 2 , X P B 1 , X R B 2 
C , X PCI) 

XNPTI,DISTI ) 

. X I,  F T I 
X R C I ) 

I/ANG/Ffi (1,1) 


/ o / 

>0) 

rd) 

F , X R 8 1 , X R B 2 
,XRB1,XRB2,XK,L0C) 


:1) 

0) 

J) 


, X R B 1 , X R B 2 


XRFDF,3) 

i 3 R , 0 ) 
IRBF, 0 ) 
IRSF, 0 ) 


) 

) 

) 

1,XRB1,XRB2,XR'I,XR2 
, X i\  C 1 ) 

, XHB  1 , XEB2 , XK , KTRIG , lEHR  ) 


G,XRCI,XK 

Xi-;FTI,DISTI  ) 
PTI,DISTI,XN 1I,DN 11) 
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c 


c 
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{ 1 , 1 ) ICODI,  XNPTi , XM 11 
>;  .ii{3)=x‘f,  'i(3)+zs(3) 
hi :■  IFF- h iO I c; H T ( x N’  1 i , x p c i , x p b 1 ^ x r 5 2 ) 

X;  (HD1?P.G:-.-2.3)G0  IQ  i75 
C.-.LL  RSPLC-  { XK  , XR2  ) 

CALL  3SCTN  ( XP,  1 , X?.2  , XRCI  , XRB  1 , XRB2  , XRD  J , XRB4  ) 
ICGDE=23 

/ I"  ITS  ( 1 , 10  : - ) ICODB  , XF.  , XR2  p XR33  . XRB  4 
MFIFr =3 
DO  95  IQ=^1,NQ 
D:-.i.R1  (10  )=-■'.  2 

D£-LP2(I0)=-3.2 
95  CCLTINUB 

D I F F R A C X 1 0 L-  3 F F 0 F.  F P S F L 3 C T 1 0 M 
T'-^w  XRFiDf  . ,c.O  , : ) GO  TO  5 
C;  .,L  SEARCH  ( 1 , 6HR5FDFHp  iRf  OF,  0 ) 

DC  110  KRr  DF=1  , N’RFDF 
CALL  READ ( TRFDF, KBAR2 , 1^3) 

C^iLL  READ  ( IPFDF  . XDL  1 p 6 , 3 ) 

CALL  READ ( lEFDF . XD32 P 6 , 0 ) 

I C 0 D S = 2 4 

WRITE{  I,  13.n  )IC0DE,KBAR2pXDB1..XDB2 
I F { XB  A R 2 . EO . K B A R 1 ) G 0 TO  1 1 0 

C ■-  L L 5 L 0 C N ( X ?.  B 3 , X F 0 4 , X F C I , X D 3 1 p X D B 2 p X K p L 0 C ) 

( LOG  . iiE~  :‘ ) GO  TO  110 

CALL  ELDPX (XR 1 p XR2p  XRCIp  XD3 1 , XDB2p  XK, RTRIG  p lERR ) 
( lERR . SO . 4 ) RETURN 
i C u u jj  — 2 u 

■■-LITE(  Ip  100  DECODE, KTRIGpXK 
I E { KTRIG  . IIQ  . 3 ) GO  TO  i ■ 3 

CALL  MR  1 ( XR  D XKpXr Cl P XNPTI . D I S TI , X N 1 1 , D N 1 1 ) 
1C0DE=26 

OPITE(1,100?)ICODE,XH  I 
ZM  10  = xM 11(3) 

X:.  1I(3)=ZN-'^  + ZS(  1 ) 

HDI?F  = H£IGHT ( XN II p XRCI , XDB 1 , XDB2 ) 

I?(HDIFf  ,GT,23. E)GO  TO  110 

CALL  SECTN(XR1 , XK , XRCI , XDB 1 , X D 3 2 , X D B 3 , X D B 4 ) 

I C 0 D E = 2 7 

lrtte(  1 . 10  ;")icode,xdb:dxdb4 

CALL  cIRpr  ( XD33  p XDBu  p XRCI,  XNPTJ  p DIST  J ) 

CALL  XH  1 ( XD33p  XDB4 , XKCIp  XIDTJ  pDISTJ  p X4  2 , DN2 ) 

I C 0 D E --- 1 2 7 

FRITS ( 1 p 1030)ICODEpXNPTJpXL2pDISTJpDN2 
HL'IFF  = HElGHT{XN2,X?.CIpXRB  i ,XRB2  ) 

XCOOS=227 

WRITE  ( 1 , 1'0..2)ICODEpHDlFF 
I? ( KDIF F . GT , -2 . 7 ) GO  TO  169 
D 0 1 g 5 1 1 = 1,  N Q 

xO=AQ+ 1 -xl 

Xri  11  ( 3 ) =ZN  : g + ZS  (IQ  ) 

H';JFF  = HEIGHI  ( XN  11,  XRCI.  XDB  DXDS2  ) 

nF.A?  = D3L  { XlDTp  XRClp  XD3  1,  XDB2p  HDIFF  p DN  1 1 ) 
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DIB'PF: 
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0DK=327 

ITF(1,10.'2)ICOD£.D£LP 
(DELF.GL.  12.5-ANr>.IQ,. 
( D E L ? , G T , D E L ? t ( I Q ) ) G 0 
( 10  . EQ  . 1 ) GO'  TO  1 1 3 
TO  105 

IQ ) =OELP 


;q.2)go 

TO  103 
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GDE 

29 

y { 

C T 1 0 


LL 


'’'It! 
I I D 


CA 
IC 
aR 
I i' 
CA 
DO 
CA 


LL  I 
CDE  = 
IXE( 

( KRF 
L ij  6 
1 U0 
L L HI 
LL  HI 


1.130 
N AFT 
0AGE  ( 
i'j  A G E ( 
I'i  A G S { 
3 3 

1,13. 
DF  . L'Q 
EAKCSt 
KHFD 
SAD  ( I 
AD  {X 


4*V  *-»  ^ -i.  . 


?)ICODE, XH2 
SK  REFLECTION 
X R 3 1 X R B 2 , X N P T X , X n p T J ) 
XRB  1 , XH32, XR  1 , XL  1 1 ) 

XRB1  ,XHd2,XR2,XR2D 

DICODS, XR;I,XK  21, XNPTI 
, 3 ) G 0 TO  14  5 
( 1 , 6HPBEDFR, IRFDF, 0 ) 

F= 1 , NRFDF 
RFDF, KBAR2, 1,0) 
RFDF,XDB1 ,6,«) 

T?  V T\  f>  O /,  ^ 


125 


CA 

IF 

If 

Cn 

X'K 
DS 
CA 
IF 
X . 

C. 

C, 


0DL=3 1 

ITE(1,13.  1)IC0DS.KBAR2,X'DB1,XDB2 

{ KBAR  2 . FQ . KBAR 1 ) GO  TO  143 

^L.  INTCPr(x?.5DXRB2,XRC.XH2I,XRB4) 

LL  3L0CN  {XF;B3,XRB4,XSC,XDD1  ,XDB2,  XJ.LOC) 
( LOG  . BQ  . -J  ) GO  TO  140 
(L0C,EQ.3)G0  TO  123 
LL  INTCPT(XR1I,XR2I,XHC,XJ,XKI) 
v(3)=ZC0R(XR1I,XR2I.XKI) 

LTA=-0 ,5 

LL  HOVE (XKI, XXI, XE  I, DELTA, I5RR ) 
(IBRR.E0.4)  RETURN 


. N E , ) G O 


TO  135 


( L 0 
TO  125 

.LL  R£PLCE(XR2I,XKI) 

,LL  INTCPT (XR1I,XPC,XDB 


I , A 


D52, XDB3 ) 


J-  - u . 


:=^32 


1 


93 


0 


c 

L 

L 


I 

a 

1 

m 

I 

■ 


"•‘,6 


L. 

c 
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W r,  I T E ( 1 , 1 E-  M 1 ) IC  0 D 2 . L C C , X J , X K I , X D 3 3 
X I'  5 3 ( 3 } “ Z C J ?.  ( X D 3 1 / X D 5 2 , X D 3 ,3  ) 

H r I r F = :•:  £ I G H I ( X R •'  I , X T B 3 . X B 3 1 # X R 3 2 ) 

Ic  (HDIFr.Gx.-2.  ’)G0  TO  i 6 'i 

C,^  LL  NR  1 ( Xr,  11 , XKI  , XRC  , XNPIJ  , DISTI  , X3  1 1 , DN  1 1 ) 

7 V r = 


X>' IK  3 ) 


NQ 


DO  130  11=1, 

I0=NQ+ 1 -II 

X;-i  11(3)  =ZX  ■ ? + Z3  ( 10  ) 

HD:.FF  = H]ijlGKT  ( XN  " I,  XPC,  XDE  1 , XDB2  ) 

Ic  (HDIFF.G-:,2:‘',".  A::0-IQ.2si,2)G0  TO  K /J 
D£LP  = D2L(XJ  iI,XRC,XDB  , X D 5 2 , HDI F F , D N 1i) 
IF ( MDIFF , £0 . 1 . AND . DEL? . GT , -3 , 2 ) GO  TO  16U 
IF f DEL? . L£ . DFLP2 ( 10 ) ) GO  TO  140 
F^D?2  ( IQ  ).=  DFLP 
CONTINUE 


6 4 


IC0DE=33 

WRITE (1,103 

") 

ICOD 

B . DELP 

IF ( DELP2 ( 2) 

.G 

E . 2 

,5)G0 

135 

CALL  REPLCE 

( XKI,  X 

E2I) 

149 

CONTINUE 

CALL  SEARCH 

( 4 

i.  6 H R 

EFDFR, 

14  5 

CALL  REPLCE 

(X 

F.2I, 

XKI  ) 

IDvR  = IT)XR3-  1 

ICODE=34 

K R I T E ( 1 , 1 0 3 

1 ) 

ICOD 

B,IDXR 

IF ( IDXR , LT . 

1 1 

) GO 

TO  ; 5 0 

f?) 


L R K 


15? 

P.  ^ 

U 

155 

IQ  = 

1 , HQ 

^ I.' 

LF ( 10, IDXR ) =AMA 

XI  (DELP  1 (IQ) 

, DELP 2 

155 

CO 

NTIN 

n -pr- 

DO 

1 6 0 

1=1 

, 3 

XI 

IG(I 

, IDXR ) =XRCI 

(1  ) 

160 

CO 

X I 

US 

GO 

TO 

165 

154 

CA 

LL  S 

i-  A ^ c 

H ( 4, 6HF. 

EFDFR , 

I R F D F , 

0) 

155 

C A 

LL  I 

i":  AGE 

( XR3 1 , X 

R 3 2 , X K I , X X ) 

159 

CA 

LL  S 

SARC 

H ( 4 , 6 H R 

X.  F i-'  r R , 

IRFDF, 

0) 

173 

CA 

LL  R 

E ? L C 

I:  ( X K , X R 

2) 

IC 

ODE=: 

3 5 

17: 


, K 


u , 


VRITE{  1 , 10  •10)ICODE,XR2 

CONTINUE 

IC0DE=36 

V.-RITS(  1 , 103e)IcODS,XR':  ,Xh,d 
;N  B’IP.RISR  factor  COFiPUTAIION 


^79 

15? 


190 


CALL  SEARCH ( 4 , 6HRSF 

HI;iG=:IDXR 

AMG=ANGLS ( VR1 , XR2 , XPC ) 
AOSr=ANG/DIST 
If { KPOS ( XNPT, XR2, Xp 1 ) . 
A ii  G 2 = A N G 1 - A N G 

cn  'pQ  -95 

ANG2=ANG 1+ANG 


C,1REF,0) 


,Q.  1 )G0  TO  190 


94 


c 


C 'J  T I-  L b 'J  T 1 0 f!  F ?.  0 H DIRECT  R A 
i ^ 5 L'  C 2 ! J 5 I Q = . , R 0 

I^{:;03(is')  .tq.:MCO  to 

r ' ?.?  = DCLF'/)  ( IQ  ) 

D 0 2 R.  0 :<  F = : , N F 

FT ( Kr, ig ) =BARFAC ( KF  # DFLP ) 

20'v  CJi;TI:njF 
ICCDF=37 

WRITE { 1 , 10  2) ICODZ , FB ( 1 . I y ) , F B ( 9 , I Q ) 
2-5  COsTINUS 

COji'T- IhUTION  FROM  REFLECTIONS 
I?  ('JIl'lG  .F.Q  ."  ) GO  TO  230 
DO  2 2 5 K I F 0 = 1 , N I M G 
f:  j 210  1=1,3 
X ( I ) =;>:i:iG  ( i , ki«g  ) 

210  C u T I :i  U E 

AhGI  = ANGLS  V XR 1 , XR2 , XRCI ) 

CALL  N3PT(XR1 ,XR2,XRCI,XN?TI,DISTI ) 

R ATI0=  ( x\NGI/DISTI  ) /ADST 
DC  220  10=1, NQ 


IF( 

nq3(iq; 

( .IQ.3 

) GO 

TO 

223 

DEL 

?=DELH i 

:iQ.Ki 

fiG  ) 

DO 

215  KF^ 

= ^,NF 

FB 

X F , I Q }- 

=F3( KF 

,IQ) 

-t-  3 A 

RFAC 

CCA 

TINUE 

223  COETINUS 


C 


T r'*  p\  "O  — . O 
J.  s.'  1 ^ 


I , F3(  1 . 1 ) ,F3(  1 , 2)  , 9;  1) 

225  CONTINUE 

COMPUTE  MFAN  ENERGY  LEVEL 

230  CALL  UR1 (XR1,XR2,XRC,XU?T,yiST,XN1,DN  I ) 

DC  292  lF=i.NF 
IP  = IF 

IF ( IF .EQ . 1 ) I?  = 5 
A = 1 3 . ^ - 1 , E - 8 4 , + =*'  I P D N 1 ) 

T 1=0  , 

DO  243  IQ=’,’!Q 

IF f NOS { IQ ) . EQ . 3 ) GO  TO  240 

NOO=NQS(IQ) 

T2  = 0 . 


F B f 9 , 2 ) 


DO  235  1=1, NQQ 

T2=T2+V3X?H ( I . IQ ) ^XLREF (I , IF, IQ ) 


235 

COF 

TI 

NU 

II 

•T'  'I 

i.  1 

+ F 

B(IF,IQ) *CQ( IQ) *T2 

240 

COR 

TI 

NU 

f-i 

Jlf 

T Tn 

A r, 

•(I 

F ) 

= XL£.  ( IF  ) + ADST*T  1 * A-FG  ( IP  ) 

242 

CON 

TI 

NU 

iT' 

il) 

ICO 

DE 

= 39 

WEITE(1,10C2)ICODE,XLE(1),XLE(9),CAP2 
A H G 1 = A 3 S { A :i  G 2 ) 
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□ 

t: 

c 

t: 


■ V) 


I 


J 

i 

1 

r 

L 

I 

E 

I 

[ 


u ^ 


•^.... 


iNl  rri 


IC0D:^  = U3 

V, -,lTi(  .'•’)ICODE,XB2,Xr<2G,XR2D,XR^0 

Ib'  { D3QR  ( XR:  , XR2  ’ ) . LT,  . . ) F.2TURN 

D2LTA-1. 

CALL  :-I0VE(X?,2,XH1  »XR  0 , DELTA  , IL'RR  ) 

I E ( I E R R , E Q . u ) RE  I U 7 L- 
IC0DS  = 4 1 

VRJTE{  1,10  ^)IC0D2,XR  . 

IF  ( D3QR  r Xin,  XR2G  ) . LT  . ■ . A- ) GO  TO  247 

Call  r s ? l c l ( x r 2 g , x r 2 ) 

0 u TO  9 0 

U7  IF (DSQR(XR7,XA2D)  ,LT.  i . l)GO  TO  25 
Call  RZPLC V ( XL  2D , XR2 ) 

GO  TO  72 

0T  f-itv4^7(5H2)CADLOI3,6F9,2/6F9.2) 

01  ?LFIlAT(5HECCDLSl3,i4,bF9.2/6F9.2) 

02  F Ar.1AT  ( 5H0OOD3(ai3, 6-  2.2) 


SL'HroijTixs  iNroR(ig) 

DIOTOJION  TL(2,9,2) 
co;-.:iON/3L:<2/NO 
co;:mon/driv  vxhph  ( 5 ) 

cn, ,nON/DRiy2/:iQS  ( 2 ) , LF 

co. :;:ou/iNXLR  •i/xlrsf(5^9^2) 

DAT:.  rL(1),  -'',(2),  TL(3),  TL(4), 


1 


j. 

r 

-•  1 

'•  ) , 

- ( 

0 \ 

\ u ; f 

T' 

L( 

25  ) , 

■:l( 

25  ) , 

T 

L( 33  ) , 

IL  { 

3u)  , 

A 

1 . 

, 7 

5 . , 3 3 , . 

48  . , 

2 . 

.57.  , 

’'87. 

t 2*  0 

2^73 

I. 

0 0 

— ■) 

05 

(IQ) 

C 

:)  V 

A ij 

' ' -1  * 

V.:t  t 

{7-: 

./30 

D 

j 

13 

T 

F=1  , 

Lt  r 

T 

I . 

? = 

T L 

(2,1 

F , I 

Q)  -T 

D 

Q 

10 

X=  1 r. 

UQQ 

X 

T '**' 

1 = 

TL(  . 

.IF 

.XQ) 

Y 

I'* 

-J  r 

(I 

, I r , 

IQ) 

= 10. 

TL(5),  TL(6),  TL(7),  TL ( 8 ) , 


1 ) , T L ( 3 2 ) , 


3 / . ^97  • ,0^.  j>55.,6o.  ,58®,  /“-O*  ,54®  ,72®  ,49. 


T P 


IQ) 


;0  CGXTIliUE 
P E T U R ,N 
7.UT) 


FOKCTION  Ai'lAG(XlV,X2V  ) 
IKD  llTGNIT'JDS  OF  VECTOR 

L'I;  E:J.3I0N'  X 1-V  ( 3 ) , X2V  ( 3 ) 
AMLG^SQRT ( DSQR ( X 1 V , X2V  ) ) 
RETURN 
END 


? 6 3 , 


3 


J 

1 


H H n n c’  LT)  tn  *jn  a'  w ro  *Tf  w w 'X  w rn  n n n i-i  :t«  m m h h *ti  w *rj  3.<  o 

nj  y,  ►.}  o o t-(  rr'  cr  r:  lo  r<<  > i?i  t"!  rt--  .v-  >■!  »-n  n ►.«  hj  itj  (;■  cj  •.:  m 


FLii'ICTION  Al’GLE  ( X 1 V,  X2V  , X3  V ) 

DIEZNSIOK  X1V(3)  ,X2V(5) ,X3V(3)  ^ 

D 3 = r SQ3  ( X V , X3V  ) 
r'23  = DS0H  ( X2V,  X3V  ) 

,Mi  G L F = 1 . 5 7 ' 8 ' - 

1F(D13^D23.FQ,: . )B£IURM 
= DSQR(X  :V,X2V) 

LE  = AC05.(D23  + D13-D'i2)/(SGRI{E13*D2J)’*'2.))  - 

U F.  N 


CTION  B.-'-.RFAC  { KF,  DEL?  ) 

FRIER  F: :T0R 
DELP.EQ.-3.2)GO  10  3 
DELP.GS.  '.2. 5)  GO  TO  4 
.<F 

KF.E0.1)IF  = 5 
ELF‘  (2.  **I?)/5.7 
.GS.78.5)G0T04 
GT  . 0 ..  AND  . A.  LT  . 7 8 . b ) GO  TO  5 
:..IQ.0.)GO  TO  6 

7.  GT,- 1.25. AND. A. LT. 0.  ) GO  TO  7 ^ 

FA  0=1.3 

URN 

iAC=4.0E-3 

URN 

FAC= (TANH ( SORT ( A ) ) *»2 ) /A/3 . 1 6 

UFN  ^ 

FAC=.315 

UR  M 

.1.  5 3 ( r.  ) ^ 

FAC=  ;TAN(SQRT{A1  ) ) ) **2/A1/5.  16 

URN 


ROUTINE  HLCCN ( X 1 V , X 2 V , X 3 V , X 4 V , X 6 V , X 6 V , LOC ) 

DATIVE  L C C A T 1 0 A OF  3 A K R I S R 

ENSIGN  X IV { 2 ) , X2V (2) , X3V ( 2 ) , X4V ( 2 ) , X5V ( 2 ) , X6V ( 2 ) , XAV ( 2 ) , XBV ( 2 ) 
= KCUT(X  ,V,X3V,X4V,X5 V) 

= I0C  + K0UT  (X2V,X3V,XUV,X5V)  ■^2 
LOC . DO , 3 ) RETURN 

L TR.KX  V.X2V,X5V,X4V,XAV,KTRI) 

LOC. EQ.-O  GO  TO  4 
K T R I . E Q . 1 ) G 0 TO  6 

2 CALL  INTC?_ ( X ^ V, X2V , X3V, X5V, XBV ) 

IF ( LOC . EG . 4 ) GO  TO  5 

3 X6V ( 1 )=XDV ( I ) 

XnV  ( 2 1 =XBV  ;■  2 ) 

RETURN 
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U I?  ( KTRI  . EQ  . ^ ) F::rUK:^ 

LCC  = U 
GO  TO  2 

5 IE ( KFQS ( X V;  , XAV , XBV ). EQ .1) GO  TO  3 

6 X dV ( 1 ) =XAV ( 1 ) 

> 5V  ( 2 ) = aAV  ■:  2 ) 

RETURN 

END 


S 1 

I UNCTION  DDL ( X 1 V , X2V, X3 V , X4V, HDIFF , Din  ) 

C FIND  PATH  LSNG,H  DIFFERENCE 

DIMENSION  >n  V ( 3 ) , X2V ( 3 ) / X3V ( 3 ) , X4V ( 3 ) ,XAV ( 3 ) , XEV ( 3 ) 
CALL  wP.FT  ( X3V  , XUV  , X ' V , XAV,  DISVA  ) 

CAIL  NRPI ( X3V, XUV, X2V, XbV, DISIB ) 

D I S TC  = D S Q R ( X R V , X A V ) 

rEL  = 3Q?.I  ( ( LISTA  + DISTB  ) * ^2  + DISTC  ) - DN  1 

II  ( H D I F F . G T , 2 . ) D E L = - D E L 
BET  U F,  N 

END 
i 1 

FUNCTION  DSOB ( X ' V, X2V ) 

DIMENSION  X 1 V ( 3 ) (.  >:2V  ( 3 ) 

DSOI  . vi 
DO  F 1=1,3 

1F  DSr.E  = ( X ^ V (I  ) - X2V  ( I ) ) * '‘2  + DSOR 

p p ■"  "I  P ?; 
r.  !.•  i V.;  i;  t'< 

L.-  At  ij 


? 1 


SUBLOUTi: 


n T ( X 1 V , X 2 V , X 5 V , X 4 V , X 5 V , X 6 V , K I R I G , I £ R H ) 


C FIND 


NEW  ENDPOi; 


DIMENSION  X n;  ( 3 ) , X2V  ( 3 ) , X3V  ( 3 ) , X 4 V ( 3 ) , X D V { 3 ) , X6V  ( 3 ) , XDUI1  ( 3 ) 
lEFI  =/! 

KTRIG=Z 
I X P I G = 1 

CALL  ?.EPLCM(X1V,XDUh) 

1 CALL  BLOCN  v X D U M , X 2 , X 3 V , X 4 V , X 5 V , X 6 V , L 0 C ) 

IF(LOC.EO,)  RETURN 

IF ( LOC . NE . 3 ) GO  TO  2 
CALL  E E P L C E ( X 2 V , X 5 V ) 

G 0 T Q 4 

2 X6V ( 3 ) =ZCCA ( X 1 V , X2V , X6 V ) 

I f { I T P. I G . E 0 . M ) GO  TO  5 
IF(AI')AG(X1V,X6V)  ,GT.C.  51  ) GO  TO  5 
ITF IG=0 

D E L T A = - 0 . 5 ' 

C ;■  L L X 0 V E ( X E U M , X D U M , X 2 V , D E L T A , I E R R ) 

GO  TO  1 

5 CELT A =-3.5 

IF( LOC.EQ. : ) GO  TO  3 

CALL  MOVE(X6V,X2V, XIV, DELTA, ISRR ) 

RETURN 


f 

•vj 
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c: 


cr 


f 1 


CALL  X0VE(Y^v,X6V,X':V,Cr:MA,IHK^) 
?:  7 r,  I G = 1 
R £ T U I-  N 
Ei'7 

FUK'CTION  H:IGri?(XlV,X2Vf  X3V,X4V  ) 


C FIND  HFIGKT  DIFFERENCE 

DIN -N SION  X1V(3),X2V(3)rX3V(3),X4V(3).XI(3) 


T I 
C 


>>  I 


CA 

LL 

1 

NICF 

I(X 

•V, 

HE 

IG 

U ^ 

— 

— i. 

0 

I ( V 

1V, 

RE 

TU 

y. 

r.  iV 

Eli 

D 

El 

c 

TION 

I 

EPS 

(X 

CH 

EC 

'< 

ON 

p 

/■  I H 

LI 

DI 

• 

KS 

X 0 N 

X 1 V ( 3 ) 

T ~ 

a. 

P S 

= 

2-' 

Dl 

ST 

z: 

A 

G (XIV 

HD 

IF 

A. 

— H 

p 

I'GH 

T ( 

D E 

12 

= 

D E 

1 

(X 

'I  V, 

LE 

L*'i 

(DE 

L 

; + D 

EL2 

IF 

( ( 

A .3  S ( D 

L2- 

DEL 

KE 

"T 

RN 

EN 

So 

r. 

0 '• 

TIN 

If'i 

A G E 

r» 

U 

j.  i"; 

AG 

E ?0 

r“  » » m 

— ;•>  j. 

DI 

N S 

ION 

X 1 V 

(3) 

AX 

=/ 

2V 

( ■'  ) 

- 

X 1 V 

{ 1 ) 

AY 

= x 

^ ; 

{ 2 ) 

- 

Y V 

( 2 ) 

AX 

V — 

A X 

aw  ^ 4- 

1 * 

* 2 

RA 

r 

0 = 

w , 

IF 

( A 

XY 

. E Q 

9 

. . ) GO 

R .A 

T1 

c = 

( { ( 

X 

3V  ( 

2 ) - 

X- 

V ( 

1 ) 

= 3 

V 

( ' ) 

+ h Y 

V . 
y. 

V ( 

2 1 

= X3V 

( 2 ) 

-AX 

x-^ 

V ( 

3 ) 

= X3 

V 

(3) 

w 

TE 

t>  ;T 

i.  ' 4.1 

D 

S. 

T-  ' 

C[f 

TIN 

E 

T 

TCP 

"T 

V • 

A-  J- 

E R C E ? 

X 

■:f 

f » 

^ vV 

D Z 

■ ? ^ 

N S 

ION 

X 1 V { 2 ) 

AX 

_ •'/ 

— A 

2V 

{ 1 ) 

- 

V 1 V ( 1 ) 

AY 

•%/ 

— A 

2V 

(2) 

- 

X ' V 

( 2 ) 

BX 

-A 

a 7 

( i ) 

- 

X5V 

( 1 ) 

J. 

— V 

— A 

U V 

( 2 ) 

- 

X 3 V 

(2) 

T 'T  ' - r»  -7-  • • f 1-  \»  ■ 

i . X A j:,  o XX  f X rt  T: . 


C1  = . Y * X 2 V ( ) - - X * X 2 V ( 2 ) 

CN-SY  *X4V ( ) _r/^x4V ( 2 ) 

I - . * X » 3 i - A i ♦ B X 
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I 


2 


S 1 
C ? 


I 


IF(D**2.LT.'.:.-6)GC  10  2 
X 5 V ( 1 ) = ( rt  X » C 2 - 3 X * C "i  ) / D 
X 5 V ( 2 ) = ( A Y * C 2 - 3 y C ',  ) / 0 
P 3 Y U >! 

D = SQRY  ( AX**2+AY-^*2  ) 

X 5 V ( 1 ) = X 1 V ( -^  ) + ( A X / D ) ♦ .3+15 

X5V ( 2 ) =X 1 V • 2 ) + ( AX/D )*  .3+15 
PEYUPN 
END 

FUNCTION  KCUT ( X 1 V, X2V , X3V , XUV ) 

Z T E P 3 T N E IF  T vv  C LINE  S E G 3 F N T 5 CROSS 

Di:. ENSIGN  X 1V  ( 2 ) , X'/V  ( 2 ) , X3V  ( 2 ) , X4V  ( 2 ) , X5V  ( 2 ) 

K C U I = / 

CALL  INTCPT ( X 1 V , X2V , X3V , X4V , X5V ) 

I F { K P 0 3 ( X 1 V , X 2 V , X 5 V ) . N E . 1 ) R E T U f N 
IF  r KF03 ( X3V , X4V, X5V ) . EG . 1 ) KCUX=  I 
PFIL’PN 
END 

FUNCTION  KPOS ( X 1 V , X2V , X3V ) 

IKD  POSITION  OF  POINT  CH  LINE 

D I E N SION  X 1 V ( 2 ) , X 2 V ( 2 ) » X 3 V ( 2 ) 

K ? 0 S = 1 

IF  ( ( (X3V  ( 1 ) - X ■’ V ( 1 ) ) -*^  ( X3V  ( 1 ) -X2V  ( 1 ) ) + ( X3V  ( 2 ) - X IV  ( 2)  ) * ( X3V  ( 2 ) - X2V  ( 2 ) 
1 ) ; ,GT.0. )KFCS=- 
P E ^ U ? M 
E V 0 

SUSrOUTINE  I'ilD?  (XlV,  X2V,  X3V) 

F I \ D C E N T E F ? 0 1 N T 

Di;.ENSICN  X^V(3)^  X2V(3),  X3V(3) 

DO  1 ?.  1=1,3 

10  X3V(I)  = ( X 1 V ( I ) +X2V ( I ) ) /2  . 


5 1 
C r 


END 

SUFTCUTINE  NIPT (XIV, X2V , X3V , X4V , DIST ) 

IND  NZAPZST  POINT  TO  LINE 

DIN  ENSIGN  X 1 V r 3 ) , X2V ( 3 ) » X3V ( 3 ) , X4V ( 3 ) , AV ( 3 ) , 3V ( 3 ) 

equivalence  { AV ( 1 ) , AX  ) , { AV ( 2 ) , a Y ) , ( AV I 3 ) , AZ ) , ( BV { 1 ) , BX ) , ( BV ( 2 ) , 3Y 

1 ) , I 5V ( 3 ) . BZ  ) 

DC  5 1=1,3 

AV  ( 1 ) =X2V ( I ) -X 1 V ( I ) 

B V f I ) = X 3 V { I ) - X 1 V ( I 1 
5 CONTINUE 

fa; 10=0. 

te;:p  = dsqp  ( :<2v , X ^ y ) 

IF ( TEMP . NE . ^ . 1 EATIC= ( A X * 3 X + A Y ^ B I + A Z * B Z ) / TEMP 
DC  10  1=1,3 

X + V (I  ) =X  1 V ( I ) +P.ATI0»  AV  ( I ) 


lOQ 


13  COi:TII!UE 

r 1 3 T = A A 3 { X L V , X 3 V ) 


IF(DI5T.SQ.^. )DI5T= 


-V^ 


R E 3 U ? X 


f 1 


S U ? r.  0 U T I N E X 5,  1 ( X 1 V , X 2 V , X 3 V , X 4 V , D I S T , X b V , D in  ) 
C FIND  NEAREST  POINT  TO  LINE  SEGilENT 


DIMENSION  X'^V(3)  ,X2V(3)  .X3V(3)  ,XUV  (3)  ,XbV(3) 
Ir (KP0S(X4V,X2V,X1V) .EO. 1 )GO  TO  2 
IF ( KP05 { X 1 V , XU  V , X2V ) . BQ . 1 ) GO  TO  4 
CALL  F.  B P L C E ( X 4 V , X b V ) 

D N 1 = D I S T 
R E E U N 

2 CALL  PEPLCE ( X 1 V, X5V ) 

GO  TO  6 

4 CALL  HBPLCE ( X2V , X5V ) 

6 DN  1=A;1AG  ( X5V,  X3V  ) 

p ■ 


I j.  I.'  JC^ 


END 


SUFxROUTINE  PEPLCE  ( X 1 V , X2V  ) 
D IE  E N S 1 0 N X 1 V { 3 ) , X 2 V ( 3 ) 

X 2V ( i } =X  1 V { ^ ) 

X2V(2)=X1V (2) 

X2V ( 3 ) =X  1 V ( 3 ) 
p E E u P :: 

end 


? 1 


S U M C u T I 4 R S E C ? N ( X 1 V , X 2 , X 3 V , X 4 V , X 5 V , X 6 V , X 7 V ) 

FIND  EFFECTIVE  BARRIE F SECTION 

E'l.lEl’SIQN  X V ( 3 ) , X2V  ( 3 ) . XjV  ( 3 i , X4V  ( 3 ) , XbV  ( 3 ) , X6V  ( 3 ) , X7V  ( 3 ) 
CALL  I N T C ? E ( X ' V , X 3 V , X 4 V , X b V , X 6 V ) 

X 6 V ( 3 ' = Z C 0 ?.  { X 4 V , X 5 V , X 6 V ) 

CALL  I 4 1 c ? ; ( X 2 V , X 3 V , X 4 V , X 5 V , X 7 V ) 

X 7 V ( 3 ) = Z C 0 ?.  ( X 4 V , X 5 V , X 7 V ) 

R '£■  '■  U ^ 

F N 


S l J B E C ' T T I N L T L I ( X 1 V , X 2 V , X 3 V , X 4 V , X 5 V , K T K I ) 

END  IF  POINT  Ii:  TRIaNGLjl  AND  LOCATE  INTERCEPT 

dx;;ei:sign  x i v ( 2 ) , X2V  ( 2 ) , X3v  ( 2 ) , xuv  ( 2 ) , Xbv  ( 2 ) 

CALL  I N T C P : ( X 1 V , X 2 V , X 3 V , X 4 V , X b V ) 

K X p j.  = 

I f ( K ? C S ( X 1 V , X 2 V , X 5 V ) . E Q . 0 ) R ET  U R N 
IF { KPOS { X3V , X5V , X4  7 ) . EG • 1 ) KTRI= 1 
PEEURN 
END 


ST 


FUNCTION  zee? ( X1 V, X2V, X3V ) 


-a 

sj 
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C FIND  7 


i; 


e;, 
e;; 
F ( 
CO 


13  zee 
:et 

n 


CCC^xDIN:- TE 

ENSIGN’  yn'(3),X2V(3),X3V(3)- 

1 = X 2 V ( M - X 1 V ( 1 ) 

2 = X2  v ( 2 ) -X V ( 2 ) 

3=X2V ( 3 ) -X1 V ( 3 ) 

IBS  ( TEMn  . GT  . ABS  ( TEX2  ) ) GO  TO  10 
?.  = X 1 V ( 3 ) + ( X3V  ( 2 ) -X  V ( 2 ) ) *TEM3/TEIi2 
U P N 

E = X 1 V { 3 ) + ( X3V { 1 ) -X  V { i ) ) «TEM3/TEM 1 


r; 


URN 


tj 


S 1 

c 


SUE ECU tine  ^OVE { X IV , X2V , X3V, deli  a , IERR  ) 
MOVE  ENDPOINT  OF  TOAD 

PIXENSION  X 1 V ( 3 ) , X2V ( 3 ) , X3V ( 3 ) 
iEp.:.=0 

ie:;?=^ahag(v^v,x3v) 

IF  f TEMP . EQ. ^ . ) GO  TO  3 
FCTF  =DELTA/TEHF 
DO  2 1=^/3 

X2’X  I I ' -X  1 V ( I)  + ( X 1 V ( I ) -X3V  ( I ) ) *FCTH 
CONTINUE 

retupn 

RET  U F N 
END 


% 1 


•>1.  ^ 
3)  I 


5.0 


FU^:CTION  I;:N(X) 
ta:.  = sin  ( X )/cos  (X) 


^ p r-'  » T T-»  1 T 

JL  i,  ^ I 


i * 


END 


FUNCTION  ACCS(X) 
Ir(X.EC.0.)  GO  TO  5 


AC0S  = ATAN ( SOAT ( ' . 2 - X * X ) / X ) 


IF  ( X . LT  . 0 , ) ACCS  = 3.  4 16  + oco; 


t;  "7>  j T 
I-  X w 


hCC^—  I ,57  2-8 


P I T U ?,  M 


END 
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USER'S  GUIDE  TO  THE  COMPUTER  PROGRAM 


4.1  Format  and  Sequence  of  Data  Input 

The  computer  program  is  written  to  accept  input  data  from  a card 
reader.  Three  types  of  input  data  are  allowed: 

Integer  — a fixed-point  number  written  without  a decimal  point. 
All  Integers  must  be  right-adjusted  within  the  allotted  field. 

Real  Constant  — a floating  point  number  written  with  a decimal 
point.  The  real  constant  may  be  situated  anywhere  within  its 
allotted  columns. 

Alphanumerie  — any  combination  of  alphabetic  and  numeric  char-- 
acters.  Alphanumeric  data  may  be  situated  anywhere  within  its 
allotted  columns. 

The  first  card  read  by  the  program  is  a title  card  with  arbi- 
trary alphanumeric  information  that  may  appear  in  columns  2 
through  60.  Up  to  five  blocks  of  data  may  follow,  in  arbitrary 
sequence,  with  information  about  the  following: 

1 Program  initialization  parameters, 

2 Road  and  vehicle  parameters, 

3 Barrier  parameters , 

4 Ground  cover  parameters, 

5 Receiver  data. 

Each  data  block  must  be  preceded  by  a control  card  containing 
an  integer  betv;een  1 and  5 in  column  5-  This  integer  char- 
acterizes the  data  block  as  listed  above.  A card  with  the 
integer  6 in  column  5 indicates  the  end  of  the  input  data 
string. 
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Multiple  cases  may  be  run  back  to  back,  whereby  each  case  starts 
with  a title  card  and  terminates  with  a card  having  a 6 in 
column  5*  The  entire  situation  need  not  be  redefined  for  each 
additional  case;  rather,  it  is  sufficient  to  specify  only  those 
aspects  of  the  problem  that  differ  from  the  case  immediately 
preceding. 

The  first  of  a series  of  cases  must  specify  initialization  para- 
meters, vehicle  and  road  parameters,  and  receiver  locations. 

If  no  information  is  provided  for  barriers  for  the  first  case, 
the  computer  program  assumes  that  there  are  no  barriers,  and 
all  calculations  related  to  barrier  diffraction  and  reflection 
are  bypassed.  Ground  absorption  is  treated  in  the  same  manner. 
Once  barrier  information  has  been  provided,  however,  failure 
to  specify  barrier  data  in  subsequent  cases  of  the  same  series 
of  cases  will  not  reinitialize  the  number  of  barriers  to  zero. 
Instead,  the  barrier  data  appropriate  to  the  preceding  case 
will  be  retained. 

Specifications  for  the  five  types  of  data  blocks  accepted  by 
the  program  are  given  below. 

4.1.1  Program  initialization  parameters 

Initialization  parameters  must  appear  in  the  first  of  a series 
of  cases.  A control  card  is  provided  with  the  integer  1 in 
column  5 to  identify  the  data  block.  Six  data  cards  .must 
follow,  in  arbitrary  sequence,  to  provide  information  about: 

1 Receiver  height  adjustment, 

2 Number  of  frequency  bands , 

3 Standard  deviation  for  cars. 
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4 Standard  deviation  for  trucks, 

5 Height  adjustment  for  cars, 

6 Height  adjustment  for  trucks. 

An  initialization  parameter  is  entered  as  a real  constant  in 
columns  1 through  10  of  each  data  card.  Column  15  contains  an 
integer  between  1 and  6 to  identify  the  parameter,  as  indicated 
in  the  above  list,  and  alphanumeric  information  is  contained 
in  columns  31  through  60.  The  last  of  the  set  of  six  cards 
must  contain  the  letter  L in  column  20  to  signify  the  end  of 
the  data  block. 

The  receiver  height  adjustment  is  a height  in  feet  to  be  added 
to  the  z-coordinate  of  all  receivers  specified  in  data  block  5* 
For  normal  conditions,  it  will  be  0.0.  For  special  cases, 
such  as  the  comparison  of  levels  at  the  ground  floor  of  a 
building  with  levels  at  a higher  floor,  a height  adjustment 
of  10.0  feet  per  floor  is  recommended.  Negative  height  adjust- 
ments are  allowable.  Such-  an  adjustment  will  not  cause  the 
program  to  assume  that  the  receiver  is  shielded  by  the  ground. 

The  number  of  frequency  bands  for  which  calculations  are  to  be 
performed  is  indicated  by  a number  (entered  as  a real  constant) 
betv/een  1.0  and  9*0.  By  using  the  number  1.0,  all  attenuations 
due  to  atmospheric  absorption,  ground  absorption,  and  barrier 
diffraction  will  be  calculated  for  a frequency  of  500  Hz. 

This  frequency  is  most  typical  for  the  spectrum  of  road  traffic 
noise.  By  using  the  number  9-0,  the  attenuations  will  be  cal- 
culated in  addition  in  eight  octave  bands  with  center  frequen- 
cies ranging  from  62.5  Hz  to  8000  Hz.  The  final  results  will 
be  more  accurate  in  this  case,  but  total  computation  time  will 
be  increased.  By  using  any  number  n between  1.0  and  9-0,  the 
attenuations  will  be  calculated  for  500  Hz  and  for  n-1  octave 
band  starting  from  the  62.5  Hz  band. 
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The  standard  deviation  for  cars  and  trucks  is  the  standard 
deviation  of  reference  levels  at  a short  distance  (r^  = 50  ft) 
from  Individual  vehicles.  The  value  recommended  for  cars  is 
2.5  (dB)  and  for  trucks  3*5  (dB). 

The  height  adjustment  for  cars  should  be  0.0,  since  tire  noise 
is  assumed  to  be  the  dominant  noise  source.  Other  height  ad- 
justments for  engine  or  fan  noise  should  not  be  made. 

The  height  adjustment  for  trucks  is  typically  8.0  (ft)  for  ex- 
haust noise.  Other  height  adjustments  for  engine  or  fan  noise 
will  result  in  slightly  higher  attenuation  by  barrier  diffrac- 
tion. A height  adjustment  0.0  for  tire  noise  is  incompatible 
with  the  noise  spectrum  stored  in  the  program. 


4.1.2  Road  and  vehicle  data 

The  control  card  for  this  data  block  must  have  the  integer  2 in 
column  5 and  a second  Integer  right-adjusted  in  columns  6 through 
10  to  indicate  the  number  of  roadways. 

The  number  of  roads  to  be  specified  depends  on  the  number  of 
traffic  arteries  with  constant  hourly  traffic  flow  that  one 
wants  to  consider.  For  a receiver  located  far  from  a multi- 
lane highway  without  ramps,  consideration  of  a single  traffic 
artery  is  sufficient.  For  receiver  locations  closer  to  the 
highway,  many  traffic  arteries  in  parallel  might  be  considered 
to  account  for  individual  lanes.  Ramps  to  highways  and  inde- 
pendent highways  are  treated  as  separate  roads. 


Follov/ing  the  control  card,  a set  of  one 
required  to  provide  vehicle  data  for  the 
card  contains  a real  constant  in  columns 
cate  vehicles  per  hour,  a real  constant 


or  more  cards  is 
first  roadway.  Each 
1 through  10  to  indi- 
in columns  11  through 
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20  to  indicate  vehicle  speed  in  miles  per  hour,  and  the  Integer 
1 or  2 in  column  25  to  indicate  vehicle  type.  The  final  card 
of  this  set  has,  in  addition  to  the  above  Information,  an  L 
in  column  31- 

An  integer  1 in  column  25  stands  for  passenger  cars,  a 2 for 
trucks.  Up  to  five  cards  for  each  vehicle  type  can  be  used  in 
random  order  to  specify  hourly  traffic  flow  at  different  average 
speeds.  Unless  the  speed  distribution  is  known  from  measure- 
ments, the  "Highway  Capacity  Manual"  is  recommended  for  deter- 
mining the  average  speed  from  the  hourly  traffic  volume  and 
the  average  highway  speed. 

Following  the  vehicle  information  for  the  first  road  are  cards 
with  Cartesian  coordinates  for  the  end  points  of  straight-line 
approximations  to  the  first  road.  The  x-y  plane  is  parallel  to 
sea  level  and  the  z-coordinate  gives  the  height  above  sea  level 
or  above  any  plane  parallel  to  sea  level.  All  coordinates  are 
in  feet.  The  x,  y,  and  z-coordinates  are  entered  as  real 
constants  in  the  fields  bounded  by  columns  1 and  10,  11  and 
20,  and  21  and  30,  respectively.  The  card  with  the  last  point 
of  the  road  has  an  L in  column  31-  A roadway  may  contain  up 
to  eleven  end  points  (i.e.,  ten  sections).  A road  containing 
more  than  ten  sections  should  be  treated  as  two  or  more  roads, 
each  having  fewer  than  ten  sections. 

Vehicle  and  road  data  are  specified  for  the  remaining  roadways 
in  a similar  manner.  Data  must  be  specified  for  as  many  roads 
as  indicated  by  the  integer  contained  in  the  second  field  of 
the  control  card. 
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4.1.3  Barrier  parameters 

The  data  block  v^^ith  information  about  obstacles  (i.e.,  "barriers") 
in  the  sound  propagation  paths  is  headed  by  a card  having  a 3 
in  column  5 and  an  integer  right-adjusted  in  columns  6 through 
10  to  Indicate  the  number  of  barriers. 

The  top  contour  of  a barrier  is  approximated  by  a straight-line 
segment,  and  no  sound  is  assumed  to  penetrate  below  this  con- 
tour. A single  barrier  may  contain  up  to  ten  sections.  The 
end  points  of  these  sections  are  specified  in  the  same  format 
as  the  end  points  of  road  sections,  except  that  the  last  point 
of  a line  is  identified  by  an  A or  an  R in  column  31  (rather 
than  L as  for  roads). 


An  R in  column  31  indicates  that 
the  top  line  of  a rigid  plane  orl 
ground,  such  as  artificial  barrie 
facades  of  buildings,  rigid  walls 

An  A in  column  31  indicates  that 
the  'top  line  of  a tilted  barrier, 
material,  an  earth  berm,  a hill, 
reflects  sound  either  v;eakly  or  t 
via  a ground  reflection. 


the  preceding  points  describe 
ented  perpendicularly  to  the 
rs'without  absorbing  material 
of  a depressed  highway,  etc. 

the  preceding  points  describe 
a barrier  with  absorbing 
or  some  other  obstacle  that 
owards  the  sky,  directly  or 


The  top  line  of  a barrier  must  not  cross  a road. 


4.1.4  Ground  cover  parameters 

The  control  card  for  this  data  block  contains  a 4 in  column  5 
and  an  integer  in  columns  6 through  10  to  indicate  the  number 
of  absorptive  ground  strips. 
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The  areas  of  ground  cover  are  described  by  the  center  line  and 
the  width  of  rectangles.  The  x,  y,  and  z- coordinates  of  one 
end  point  of  the  center  line  are  given  as  real  constants  in 
the  fields  between  columns  1 and  10,  11  and  20,  and  21  and  30, 
respectively  of  a single  card.  The  same  card  contains  the 
width  of  the  rectangle  in  real  formats  in  the  field  between 
columns  31  and  ^0.  The  x,  y,  and  z- coordinates  of  the  other 
end  point  of  the  center  line  are  written  on  the  next  card 
together  v/ith  a G or  a T in  column  31  • 

A G identifies  the  ground  cover  as  high  grass  or  shrubbery, 
and  a T identifies  the  ground  cover  as  trees. 

The  rectangles  with  ground  cover  must  not  cross  a road. 

4.1.5  Receiver  data 

The  control  card  contains  a 5 in  column  5 and  an  Integer  in 
columns  6 through  10  to  indicate  the  number  of  receivers.  A 
card  is  provided  for  each  receiver  to  Indicate  the  x,  y,  and 
z-coordinates  of  the  receiver  location.  These  data  are  entered 
as  real  constants  in  the  fields  bounded  by  columns  1 and  10, 

11  and  20,  and  21  and  30,  respectively. 

A maximum  of  fifteen  receivers  is  allowed. 

A receiver  cannot  be  located  on  a road,  nor  on,  over  or  under- 
neath a top  line  of  an  obstacle  (barrier),  nor  on  a ground 
strip.  However,  a receiver  can  be  located  betv;een  two  adjacent 
ground  strips  if  the  location  is  not  Identical  with  an  end 
point  of  the  ground  strip  center  line. 
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4.2  Error  Messages 

Errors  are  detected  by  the  computer  in  the  follov/ing  cases. 

a.  If  the  top  line  of  an  obstacle  intersects  a roadv;ay,  the 
computer  prints  out 

ILLEGAL  BARRIER  INTERSECTS  ROADWAY 

together  with  the  following  data:  receiver  number,  roadway 

number,  road  section  number,  barrier  number,  barrier  section 
number . 

The  computer  then  proceeds  with  the  next  case. 

b.  If. the  center  line  of  a ground  strip  intersects  a roadway, 
the  computer  prints  out 

ILLEGAL  GROUND  STRIP  INTERSECTS  ROADWAY 

together  with  the  following  data:  receiver  number,  roadway 

number,  road  section  number,  ground  strip  number. 

The  computer  then  proceeds  with  the  next  case. 

c.  If  the  number  of  reflections  contributing  to  the  sound  level 
at  a receiver  exceeds  10,  the  computer  prints  out 

TOO  MANY  REFLECTIONS 

together  v;ith  the  following  data:  receiver  number,  roadway 

number,  road  section  number. 

The  computer  then  proceeds  with  the  next  receiver. 

d.  Should  the  geometry  of  the  situation  be  such  that  a line 
segment  that  is  needed  for  computation  appears  to  have  zero 
length,  the  program  continues  with  the  next  receiver  with- 
out giving  an  error  message. 
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4.3  Data  Output 

The  data  output  starts  immediately  after  the  input  file  is  read. 
The  heading  TRAFFIC  NOISE  PREDICTION  is  typed  and  input  data 
are  then  typed  in  the  following  order: 

Title  card; 

Program  initialization  parameters  (for  the  first  of  a series  of 
cases  only ) ; 

Vehicle  and  road  parameters,  where  all  parameters  of  Type  1 
vehicles  (passenger  cars)  on  a given  road  are  typed  first, 
following  by  parameters  for  Type  2 vehicles  (trucks); 

Barrier  parameters,  with  the  type  of  barrier,  A or  R,  given 
in  parentheses  after  the  barrier  number; 

Ground  cover  parameters,  with  identification  G or  T following 
the  ground  strip  number  in  parentheses; 

Receiver  data. 


The  title  card  is  then  printed  out  again  together  with  a- head- 
line for  the  receiver  identification.  The  output  consists  of 
the  following  data  for  each  receiver.  After  the  sound  pressure 
level  for  a given  receiver  has  been  calculated,  the  receiver 
number  and  coordinates  are  typed.  The  octave  band  center  fre- 
quencies are  typed  out;  underneath  them  are  the  calculated  A- 
weighted  octave  band  levels.  The  final  typeouts  consist  of 
the  A-weighted  overall  level  based  on  attenuations  at  500  Hz, 
which  appears  under  the  heading  A(500),  and  the  A-weighted 
overall  level  resulting  from  the  octave  band  calculations, 
v/hich  appears  under  the  heading  A. 


Ill 


4.4  Recommendations  for  Usage 


The  computation  time  required  to  compute  sound  levels  at  a given 
receiver  increases  with  the  number  of  roadways,  barriers,  and 
ground  strips,  but  it  is  independent  of  the  total  number  of 
receivers.  Therefore,  if  one  wishes  to  analyze  a situation 
in  which  many  receivers  and  barriers  are  located  along  or 
near  a highway  many  miles  in  length,  we  suggest  that  the  prob- 
lem be  modularized  so  that  a smaller  geographical  area  can  be 
associated  with  each  receiver. 


Maps  showing  the  elevation  of  the  terrain  should  be  used  for  a 
detailed  description  of  the  roadway  and  of  top  lines  of  hills 
(barriers  of  Type  A) . It  should  be  borne  in  mind  that  the 
computer  program  accounts  for  5 dB  attenuation  for  sound  rays 
grazing  over  hills. 

Results  obtained  w'hen  ground  absorption  is  considered  should  be 
interpreted  v;ith  caution.  Because  of  the  lack  of  reliable  data 
in  this  area,  results  from  the  first-order  approximations  used 
in  the  computer  program  are  intended  for  comparison  with  field 
data  rather  than  for  prediction  purposes. 

4.5  Test  Cases 

As  an  example,  a listing  of  input  data  for  two  cases  and  the 
corresponding  computer  printout  is  given  on  the  following 
pages.  A graph  of  the  geometry  in  the  x-y  plane  is  shovm  on 
Fig.  13. 

In  the  example,  test  case  1 involves  two  roadways,  a tree  zone, 
and  two  receivers.  The  curved,  major  roadway  with  1000  passen- 
ger cars  per  hour  at  an  average  speed  of  60  mph  and  with  100 
trucks  per  hour  at  50  mph  is  approximated  by  two  straight 
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FIG. 13  PLAN  VIEW  OF  THE  GEOMETRY  USED  IN  THE  TEST  CASES 
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sections.  The  straight,  minor  roadway  with  30D  passenger 
cars  per  hour  at  an  average  speed  of  ^0  mph  has  only  one 
section. 

Since  most  of  the  sound  is  radiated  from  the  major  roadway  and 
the  receiver  locations  are  symmetrical  about  its  first  section, 
the  sound  pressure-  levels  calculated  at  the  two  receivers  are 
about  equal. 

Test  case  2 involves  the  same  parameters  as  in  test  case  1 as 
well  as  tv;o  barriers.  One  of  the  barriers  might  be  a hill  or 
an  earth  berm  and  is  therefore  of  an  absorptive  type.  The 
other  barrier  might  be  a building  and  is  specified  by  four 
reflective  sections. 

The  absorptive  barrier  is  between  10  and  20  ft  higher  than  the 
major  roadway  and  shields  receiver  no.  1 against  noise  from 
almost  the  entire  major  roadv;ay.  However,  diffraction  of 
sound  over  the  barrier  and  reflections  from  the  other  barrier 
limit  the  noise  reduction  at  this  receiver  to  about  8 dB(A). 
Receiver  no.  2 is  only  partially  shielded  by  the  reflective 
barrier  against  noise  from  the  major  roadway  so  that  the  noise 
reduction  is  only  about  3 dB(A) . 
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